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Abstract  (cont.) 


The  NORSAR  Detection  Processing  system  has  operated  throughout  the 
reporting  period  with  an  average  uptime  of  98.9  per  cent.  A  total  of 
1940  seismic  events  have  been  reported  in  the  NORSAR  monthly  seismic 
bulletin.  The  performance  of  the  continuous  alarm  system  and  the  auto¬ 
matic  bulletin  transfer  by  telex  to  AFTAC  have  been  satisfactory.  Pro¬ 
cessing  of  requests  for  full  NORSAR/NORESS  data  on  magnetic  tapes  has 
progressed  according  to  established  schedules. 

A  new  Detection  Processing  system  has  been  Implemented,  providing 
Improvements  especially  in  the  data  quality  checking  and  reduction  in 
false  detections  caused  by  occasional  data  spikes.  The  NORSAR  beam 
table  and  event  procetsing  algorithms  remain  unchanged.  Processing 
code  has  been  developed  for  on-line  regional  array  analysis  based  on 
SUN  3  computer  systems. 

Field  maintenance  activity  has  Included  regular  preventive  maintenance 
at  all  subarrays  and  occasional  corrective  actions  when  required.  No 
special  problems  have  been  noted  in  the  performance  of  field  systems. 
In  addition,  field  maintenance  personnel  have  participated  in 
establishing  the  new  regional  array  in  Finnmark,  Northern  Norway. 

A  study  of  Pn  scattering  has  been  conducted,  using  as  a  data  base 
NORESS  recordings  of  six  mining  explosions  in  the  Blasjp  area,  at  an 
epicentral  distance  of  about  300  km.  Wideband  f-k  analysis  of  these 
event  recordings  show  a  pronounced  and  consistent  anomaly  in  the  fre¬ 
quency  range  2-4  Hz,  which  is  also  the  frequency  band  of  maximum  Pn 
energy.  This  is  found  to  be  consistent  with  scattering  by  topographic 
relief  of  the  Moho  discontinuity,  and  the  interred  location  of  the 
proposed  topographic  feature  coincides  with  the  border  of  the  Oslo 
Graben. 

A  study  of  the  effect  of  lateral  heterogeneities  on  short  period  sur¬ 
face  waves  (T  <  20  s),  using  numerical  modelling  with  a  coupled-mode 
scheme,  has  been  initiated.  So  far  ft  has  been  found  that  the  coupled- 
local  modes  method  can  be  applied  to  model  wave  propagation  in 
realistic  structures,  needing  a  reasonable  amount  of  lateral  zone;, 
local  modes  to  achieve  satisfactory  precision  in  transmission  and 
reflection  coefficients.  The  plan  is  now  to  apply  this  method  to 
laterally  heterogeneous  structures  around  Scandinavia,  and  compare  the 
results  of  the  numerical  modelling  to  surface  wave  data  recorded  at 
NORSAR  and  NORESS. 
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(  The  problem  of  threshold  determination  for  regional  event  detection 

has  been  investigated,  addressing  in  particular  the  tradeoff  between 
detector  sensitivity  and  false  detections  caused  by  noise  fluctuations 
or  SLOchastic  nature.  It  has  been  found  that  for  the  same  false  alarm 
rate,  thresholds  of  incoherent  beams  can  be  set  at  A. 5-7  dB  lower  than 
for  coherent  beams.  This  implies  that  for  Pn  detection,  coherent  beams 
will  usually  outperform  incoherent  beams,  since  the  SNR  gain  is  typi¬ 
cally  8-15  dB.  However,  for  secondary  phase  detection,  the  incoherent 
beams  will  usually  be  the  best,  since  the  coherent  beam  gain  in  these 
cases  is  much  lower.  A  particularly  promising  result  is  that  inco¬ 
herent  beamforming  based  on  horizontal  channels  appears  to  provide 
more  effective  noise  variance  reduction  than  is  obtained  using  ver¬ 
tical  channels.  Also  a  large  number  of  noise  detections  can  be  iden¬ 
tified  as  such  using  the  slowness  estimates  provided  by  a  small- 
aperture  array. 

Initial  data  recorded  at  the  new  regional  array  in  Finnmark,  Northern 
Norway,  have  been  processed  and  compared  to  NORESS  recordings.  The 
Finnmark  array  is  found  to  provide  the  same  excellent  (better  than  /N) 
Twice  suppression  as  NORESS  in  selected  filter  bands,  using  "optimum" 
subgeometries  for  beamforming.  The  background  noise  level  at  Finnmark 
is  even  lower  than  that  of  NORESS  at  high  frequencies  (greater  than 
2  Hz),  but  somewhat  higher  than  for  NORESS  at  the  low  frequency  end, 

'  due  co  closer  distance  from  the  coast.  A  feature  of  Finnmark  recor- 

t  dings  not  seen  at  NORESS  is  the  presence  of  strong  Rg  phases  for 

mining  explosions  out  to  about  350  km  epicentral  distance.  A  two-week 
test  period  of  joint  processing  of  the  two  arrays  has  shown  that  the 
large  majority  of  small  events  (Mg<2.0)  are  detected  by  only  one 
'  array,  whereas  larger  events  (ML>2.0)  detected  by  either  array  in  most 

cases  have  at  least  one  confirming  phase  detection  from  the  other 
array . 

Analysis  of  P  coda  and  Lg  phases  recorded  at  NORSAR  for  a  set  of  72 
large  explosions  from  the  Shagan  River  test  area  near  Semipalatinsk , 

(  USSR,  has  shown  that  such  data  can  provide  event  magnitude  estimates 

comparable  in  stability  to  world-wide  network  mb-  A  systematic  dif¬ 
ference  in  m5(P)-m{,(bg)  has  been  observed  for  the  northeast  and  south¬ 
west  portion  of  the  Shagan  River  test  area.  There  are  indications  that 
this  is  mainly  due  to  a  bias  in  m[,(P)  estimates,  and  if  so,  the  m^lg) 
values  observed  at  NORSAR  would  be  expected  to  give  more  reliable 
relative  yield  estimates  than  can  be  obtained  from  P-wave  data. 

A  large-scale  study  comprising  NORESS  recordings  of  101  events  from 
the  North  Sea/western  Norway  region  has  been  conducted  with  the  aim  to 
determine  optimum  beam  deployment  for  regional  phase  detection.  The 
study  has  reconfirmed  that  a  number  of  narrow  band  filters  applied  in 
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parallel  and  In  combination  with  array  beams  based  on  filter-dependent 
subgeometries  produce  the  best  detection  performance.  In  particular, 
the  study  has  demonstrated  the  usefulness  of  high-frequency  (10-16  Hz) 
coherent  beams  for  Pn  phase  detection.  It  has  furthermore  been  shown 
that  high  frequency  incoherent  beams  (10-16  Hz)  formed  from  horizontal 
channels  significantly  enhance  detectability  of  Sn  phases.  The  best 
results  in  Lg  phase  detection  are  found  for  low  frequency  incoherent 
beams  (1.5-4. 5  Hz)  combining  both  horizontal  and  vertical  channels. 
Further  work  will  in  particular  be  directed  toward  analyzing  events 
from  other  regions  in  order  to  determine  a  beam  deployment  that  provi¬ 
des  optimized  detection  in  the  entire  regional  distance  range. 

The  local  geology  of  the  two  regional  array  sites  in  Norway  nas  been 
studied,  with  emphasis  on  compiling  information  that  ir  considered 
relevant  from  a  seisniological  viewpoint.  At  this  stage,  the  main  pur¬ 
pose  has  been  to  present  results  of  various  surveys  conducted  by  the 
Norwegian  Geological  Survey  and  NORSAR  scientists.  In  the  future,  this 
information  is  expected  to  become  useful  in  explaining  observed 
features  in  the  recorded  seismograms  from  the  two  arrays. 
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I.  SUMMARY 

This  Semiannual  Technical  Summary  describes  the  operation,  maintenance 
and  research  activities  at  the  Norwegian  Seismic  Array  (NORSAR)  for 
the  period  1  April  -  30  September  1987. 

The  uptime  of  the  NORSAR  online  detection  processor  system  has 
averaged  98.9  per  cent  during  the  reporting  period  as  compared  to  99.6 
for  the  previous  period.  04C  and  06C  have  been  the  most  frequently 
affected  systems  due  to  cable  damage  (04C)  and  synchronization 
problems  (06C).  A  total  of  1940  seismic  events  were  reported  In  the 
NORSAR  monthly  seismic  bulletins  in  this  period,  giving  a  daily 
average  of  10.6  events. 

The  new  Detection  Processor  (DP)  was  put  into  routine  operation  on  day 
289.  The  detector  program  is  executing  on  the  IBM  4341/L01  computer 
for  the  NORSAR  array.  There  are  no  changes  in  the  beam  table.  From  day 
320,  the  MODCUMP  processor  has  been  doing  data  acquisition  only,  as 
the  beamforming  and  filtering  tasks  have  been  moved  to  the  IBM  4341. 
The  new  detection  process  has  lead  to  a  significant  reduction  in  load 
on  the  event  processor,  due  to  a  g'°atly  reduced  number  of  "spike" 
detections.  There  are  no  changes  in  the  event  processor  code.  There 
are  also  no  changes  in  the  'online'  RONAPP  detection  processing  for 
the  NORESS  array,  but  the  new  detector  package  now  being  used  for 
NORSAR  has  also  been  applied  to  NORESS  data  and  is  executing  regularly 
on  tapes  with  data  from  the  new  Finnmark  regional  array. 

A  new  program,  SEGPRO  -  Segmented  processing,  has  been  developed  to  do 
the  onset  determination  as  well  as  phase  velocity  and  azimuth  measure¬ 
ments  on  detections  reported  from  the  new  detector  package,  on  a 
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routine  basis.  A  new  event  processor  code  is  under  development  to  pro¬ 
cess  NORESS  and  Finnmark  array  data. 

The  data  acquisition  and  real-time  processing  systems  tor  the  new 
Finnmark  array  will  be  based  on  three  Sun  3  computers.  The  IBM/PC 
currently  used  for  controlling  critical  tasks  in  the  NORSAR/NORESS 
systems  will  also  be  used  for  the  new  array. 

Much  time  and  effort  has  been  spent  by  the  field  maintenance  personnel 
on  preparatory  work  in  connection  with  the  new  array  in  Finnmark,  and 
the  actual  installation  of  the  array  required  extensive  field  work 
during  the  summer  of  1987. 

The  r«f«arc-  activity  Is  summarized  in  section  Vll.  Section  VII. 1 
presents  a  scattering  model  of  regional  Pn  wave  propagation.  Section 
VII. 2  describes  preliminary  tests  for  surface  waves  in  2-D  structures. 
In  Section  VII. 3  false  alarm  statistics  and  threshold  determination 
‘oi  regional  event  detection  are  presented.  Section  VII. 4  gives  ini¬ 
tial  results  from  analysis  of  data  recorded  at  the  new  regional  array 
in  Finnmark.  In  Section  VII. 5  the  potential  of  using  NORSAR  recorded  P 
coda  and  L,g  waves  to  obtain  stable  magnitude  estimates  of  large  Semi- 
palatinsk  explosions  is  investigated.  Section  VII. 6  describes 
Investigations  of  optimum  beam  deployment  for  NORESS  using  a  North  Sea/ 
western  Norway  data  base.  In  Section  VII. 7  the  local  geology  of  the 
regional  array  sites  in  Norway  is  presented. 


II. 


OPERATION  OF  ALL  SYSTEMS 


II. 1  Detection  Processor  (DP)  Operation 

There  have  been  88  breaks  in  the  otherwise  continuous  operation  of  the 
NORSAR  online  system  within  the  current  6-month  reporting  interval. 

The  uptime  percentage  for  the  period  is  98.9  as  compared  to  99.6  for 
the  previous  period. 

Fig.  II. 1. 1  and  accompanying  Table  II. 1.1  both  show  the  daily  DP 
downtime  for  the  days  between  1  April  and  30  September  1987.  The 
monthly  recording  times  and  percentages  are  given  in  Table  II.  1.2. 

The  breaks  can  be  grouped  as  follows: 


a)  Hardware  failure  9 

b)  Stops  related  to  program  work  or  error  0 

c)  Hardware  maintenance  stops  6 

u )  Power  jumps  anti  breaks  5 

e)  TOD  error  correction  6 

f)  Communication  lines  67 


The  total  downtime  for  the  period  was  50  hours  and  26  minutes.  The 
mean-time-between-failuras  (MTBF)  was  2.5  days,  as  compared  to  2.0  for 
the  previous  period. 


J.  Torstveit 


*  Mean-t lme-between-failures  -  total  uptime/no.  of  up  Intervals 


e  II. 1.2  Online  system  performance,  l  April  -  30  September  1987 


DP  uptime 

DP  uptime 

No.  of 

No.  of  days 

DP  MTBF* 

hours 

Z 

DP  breaks 

with  breaks 

(days) 

-  6  - 


I 

w 


W  X  w  w 

X  o«« 

^33^ 

i  PSOSQ^DMMOS 

$33S5S!££3 
5£2amiSS£2 

KKKOO 
J  u  w  u  o  u  w 
^22050002 
<hhQOOOOh 

jjjpiPssd 

)OON(M^OMno«N( 


mUUU  t*  -J  CU 

h  «s  a  sc 

'323S..333 

m  hKhhhh 

<x<ea<<<< 
u.  o  u.  u,  u,  u.  u. 
a 


<  <  <  <  CL  G.  - 
kh.b.|i.SZl 

Op  ■  : 

U2WWWUOWXWWWWWW 
22Z2QQ20!ZS2Z22 
mmmmOO*- 

J  >J  >J  iJ  Z  £  k 


<OM'NH(NnvOV'OlMHOi«MnOON«rOnift'0®Vlfl<0'f 
04  in  04  n  04  H  *t  HVNOIftH  CM  PI  Ifl  ^  04  V  M  «f  PI 


^OfMN'Oin'Ocoo*ONO'0'fiOcOii«oi'Onor>'fion(M««rc 


HHHHHNOin'OOCO'rVO'O'OCODCOfflOHNrrinOCOOOOO 

ooooHHHofojoiMnw^iniiiirtin'OoooooocooacD 


a 

$ 

o 

PO 

I 


w  £ 
w  x 
u  o 
uu 
o 
fi 

04 

04  04 

DO 

H 

2  0) 


04  (0 

X  X 

O  'C 

JUtdUUUUUUUOblUUU  UUUUUUUU 

iaaaficsKKacsQSKKBiDaaKKKaaa 

‘133333333§3333S333P3333 

M  W  W  M  M  M  M  M  M  M  ►-<  »H  M  M  M  M  HHHhShHH  ~i  M  M  M 
<<<<<<<<<<<<<<<<W<<<<H<<<<<<<< 

o.u^fc.u^u.u.u.u.u.u.^u^u.ou.fr.u.u.wu.u.u.u.u.u.u.u. 

UUUUUUUUUUtdUUUUUXilUUU  UUMWMMUU 
ZZZZZZZZZZZZZZZZZZZZZQZZZZZZZZ 


2 

< 

W  H 

£3 

u.  w 
o 


5* 


U)>4  HHrtHHHNPmPt>nr«MOO'«JiO'0'0'OOOHHHHr4f4H 

N  <  (^(JlOlOO'O'l^O'lT'OO'O'C'IJ'O'ff'O'C'OCJ'OOOOOOOOOO 

JO  rlHHHHHHHrtH 


Table  IX. 1.1  Daily  DP  down  ime  in  the  per 
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Table  II. 1.1  (cent. 


II. 2 


Array  communications 


Table  VII. 2.1  reflects  the  systems  performance  in  the  reporting 
period . 

04C  and  06C  have  been  the  most  frequently  affected  systems  due  to 
cable  damage  (04C)  and  synchronization  problems  (06C). 

The  communication  machine  (Modcomp)  has  failed  two  times,  in  May  and 
September,  a  voltage  and  a  "Floating  point  processor"  failure, 
respectively.  The  latter  failure  caused  problems  during  the  restarting 
phase. 

Summary 

Apr  In  April  most  systems  were  affected  especially  in  connection 

with  our  attempts  to  localize  the  noise  source  common  to  all 
communication  channels.  Systems  were  switched  off  one  by  one, 
but  also  in  groups  of  2  and  3-  Later  NTA/Hamar  and  Lillestrpra 
did  observations  which  revealed  noise  in  the  transmission 
channel  equipment  between  H.tmar  and  Oslo.  NTA  switched  to  a 
spare  channel  group  which  for  a  period  gave  us  improved 
performance. 

May  Apart  from  06C  which  was  affected  weeks  19,  21  and  22  causing 

0.5%  outages,  the  systems  were  most  reliable  in  May.  The  high 
error  figures  were  caused  by  a  Modcomp  failure  week  19, 
representing  an  outage  of  approx.  15.4%  that  week.  Figures  in 
brackets  represent  system  performance/reliability  exclusive 
Modcomp  failure  (see  Table  II. 2.1). 
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Jun 


Jul 


Aug 


Sept 


In  spite  of  varying  line  quality  at  06C  weeks  25  and  26  and 
relatively  high  error  figure  (9.5%)  this  period  may  be 
characterized  as  satisfactory  for  the  remaining  systems, 
although  a  0.21%  error  figure  at  04C  could  have  been  better. 

Defective  cable  between  Berge  Central  and  the  04C  CTV  was  a 
problem.  Roadwork  was  going  on  in  the  04C  area,  affecting  the 
communication  cable  which  had  to  be  cut/spliced  at  different 
places  along  the  path.  NTA,  who  was  responsible  for  the  cable 
work,  assumed  (by  the  end  of  the  month)  that  it  would  take 
another  2-3  weeks  before  the  work  could  be  finished. 

02C  (8.2%  outages)  was  frequently  affected,  but  NTA  did  not 
localize  the  irregularity. 

06C  communication  system  had  synchronization  problems  which 
may  have  been  a  result  of  a  "stressed"  Modcomp,  caused  by 
04C . 

High  performance  by  the  remaining  systems. 

Week  35  NTA  finished  the  cable  work  at  04C,  and  the  error 
figures  decreased  to  0.06%  that  week. 

NTA  was  told  to  check  the  06C  communication  line  thoroughly. 
Outages  occurred  frequently,  but  the  station  was  back  in 
operation  after  restarting  the  Modcomp. 

Most  satisfactory  performance  of  the  other  systems. 

After  the  cable  work  finished  week  34,  04C  communication 
system  functioned  most  satisfactorily  the  three  following 
weeks  (35,  36,  37),  but  16  September  data  again  disappeared. 
After  a  temporary  repali  the  system  was  back  in  operation  23 
September.  Between  24  and  28  September  the  performance  was  bad. 


06C  pattern  did  not  change  compared  to  July/August. 
Reaynchronlzation  is  possible  and  brings  the  system  up  again. 

The  NMC  staff  will  as  soon  as  they  have  finished  the  installa¬ 
tion  of  the  new  array  in  Finnmark  check  relevant  equipment 
(power,  modem,  Loop  Control,  Logic,  etc.),  which  may  be 
carried  out  in  conjunction  with  NTA/Hamar  and  the  NDPC  staff. 

O.A.  Hansen 
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III.  ARRAY  PERFORMANCE 


Event  Processor  Operation 

In  Table  III-l  some  monthly  statistics  of  the  Event  Processor  opera¬ 
tion  are  given: 


Teleseismic 

Core  Phases 

Sum 

Daily 

APR  87 

264 

69 

333 

11.1 

MAY  87 

254 

62 

316 

10.2 

JUN  87 

275 

60 

335 

11.2 

JUL  87 

332 

85 

417 

13.5 

AUG  87 

205 

57 

262 

8.5 

SEP  87 

210 

67 

277 

9.2 

1540 

400 

1940 

10.6 

Table  III-l  Event  Processor  Statistics,  April  -  September  1987. 


B.  Paulsen 
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IV  IMPROVEMENTS  AND  MODIFICATIONS 

IV.  1 _ NORSAR  detection  processing 

The  'new'  NORSAR  detection  processor  was  put  into  routine  operation  on 
day  289,  after  a  long  period  of  parallel  processing  for  verification. 

The  change  to  the  new  detection  processor  was  done  under  a  period  of 
bad  data  conditions.  Three  subarrays  were  temporarily  out  of  order, 
and  many  communication  errors  were  present.  However,  these  bad  data 
conditions  were  exploited  to  gain  valuable  experience  with  the  data 
quality  control  part  of  the  system.  We  now  have  a  program  system  that 
controls  all  bad  data  conditions  due  to  communication  errors  and 
resynchronizations.  Moreover,  the  data  quality  control  procedure  is 
changed  from  treating  only  the  special  features  of  the  NORSAR  array  to 
a  more  general  parameterized  control  routine.  This  routine  performs 
data  control  on  the  NORSAR,  NORESS  and  new  Flnnmark  arrays  using  the 
same  program  code  for  all  arrays. 

The  output  from  the  detector  system  is  in  the  form  of  accumulated 
sequential  files.  One  file  gives  a  list  of  detections,  another  gives  a 
list  of  'events'  from  the  data  quality  control  part.  Data  quality 
control  reporting  Includes  time  and  length  of  data  gaps  for  individual 
channels  and  for  the  array  due  to  communication  errors  and/or 
recording  gaps.  Moreover,  spikes  are  reported  if  detected,  using  max. 
amplitude  comparison  between  one  sensor  and  the  average  or  median 
value  of  all  sensors.  (Other  parameterized  procedures  may  also  be 
selected).  Using  this  procedure,  spikes  of  significant  amplitudes  are 
effectively  removed. 

The  detector  program  is  executing  on  the  IBM  4341/L01  computer  for 
the  NORSAR  array.  There  are  no  changes  in  the  beam  table.  Coherent 
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beams  1  -  180  and  Incoherent  beams  1-64  are  directed  as  before 
(Beam  tables  451  and  271,  respectively). 

To  reduce  computer  power  requirements,  the  band  pass  filtering  is 
done  using  a  20  Hz  sample  rate,  whereas  the  array  beamforming  is  done 
using  10  Hz  sampling.  This  reduction  in  sample  rate  has  in  practice 
shown  no  degradation  of  the  signal-to-noise  ratio  within  the  frequency 
bands  used.  We  even  get  a  better  performance  relative  to  the  'old' 
system  where  the  beamforming  was  done  using  two  steps:  Firstly, 
subarray  beamforming  was  done  using  a  small  set  of  beams,  and  then 
array  beamforming  was  done  using  this  set  of  subarray  beams-  The  new 
detect or  forms  array  beams  using  all  sensors,  with  individual  delays 
for  all  beams,  thus  obtaining  a  better  signal-to-noise  ratio. 

The  4341  computer  executes  M0DC0MP  communication,  data  recording  task, 
detection  processing,  event  processing  and  AFTAC  telex  message 
handling  for  the  NORSAR  array.  For  60  seconds  of  data,  the  computer 
spends  on  the  average  50  seconds  of  elapsed  time  to  execute  all  these 
processes.  The  detection  processor  runs  in  a  mode  where  it  takes  a 
3-minute  break  every  time  it  catches  up  with  real  time. 

In  this  mode,  close  to  real  time,  the  program  may  execute  for  several 
weeks  without  interruptions.  For  this  reason,  the  detector  system  has 
a  built-in  task  which  may  communicate  externally,  and  performs  file 
management  without  interrupting  the  continuity  in  the  data  or  the  exe¬ 
cution  of  the  program.  However,  the  detector  program  has  so  far  been 
stopped  every  7  days  to  do  manual  maintenance  of  output  files  and 
backup.  (The  "old"  online  disk  loop  with  data  and  detection  lists  is 
maintained. ) 
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IV. 2 _ MODCOMP  subarray  communication 

From  day  320,  the  MODCOMP  processor  has  been  doing  data  collection 
only,  as  subarray  beamforming  and  filtering  was  removed  from  the 
code,  'the  removing  of  cubgrray  beamforming  end  filtering  has  relieved 
the  CPU  comsumption  on  the  MODCOMP,  leading  to  a  much  better  perfor¬ 
mance  on  handling  the  SLEM  communication. 

IV.  3 _ NORSAR  event  processing 

The  new  detection  process  has  lead  to  a  significant  relief  on  the 
event  processor,  as  most  of  the  'bad  data  events'  are  removed. 

There  are  no  changes  in  the  event  processor  code. 

An  actual  improvement  in  automatic  processing  results  requires  a  new 
time  delay  correction  data  base. 

IV. 4 _ NORESS  detection  processing 

There  are  no  changes  in  the  'online'  RONAPP  detection  processing,  but 
the  new  detector  package  now  being  used  for  the  NORSAR  array  pro¬ 
cessing  has  also  been  applied  to  NORESS  data  and  gives  results  inden- 
tical  to  those  of  RONAPP.  Potential  gain  in  detection  capabilities 
using  extended  beam  sets  for  NORESS  are  discussed  elsewhere  in  this 
report . 

IV. 5 _ Flnnmark  regional  array  detection  processing 

The  detector  package  now  being  used  lor  the  NORSAR  array  data  is  now 
executing  regularly  on  tapes  with  data  from  the  new  Finnmark  regional 
array.  Resu'ts  from  this  detection  processing  are  given  elsewhere  in 
this  report. 
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IV. 6 _ NORESS/Finnmark  array  event  processing 

There  are  no  changes  in  the  event  processing  for  NORESS  data.  The  event 
processing  consists  of  phase  identification,  phase  association  and 
event  epicenter  determination,  and  subsequent  plotting  of  relevant 
data  and  f-k  spectra. 

A  new  program,  SEGPRO  -  Segmented  processing,  has  been  developed  to  do 
the  onset  determination  as  well  as  phase  velocity  and  azimuth  measure¬ 
ments  on  detections  reported  from  the  new  detector  package,  on  a 
routine  basis.  The  SEGPRO  package  is  regularly  processing  Finnmark 
array  data.  The  detector  and  the  SEGPRO  package  together  perform  all 
RONAPP  functions,  and  are  using  the  same  data  Interface  routines.  This 
means  that  the  programs  may  use  data  from  disk  loops,  tape,  disk  files 
or  any  other  media  or  format  that  is  supported. 

A  new  event  processor  code  is  under  development  to  process  NORESS 
and  Finnmark  array  data  using  these  routines.  The  processor  will  use 
the  SEGPRO  output  files  for  detection  parameters  and  f-k  matrices, 
thus  preparing  for  a  complete  replacement  of  RONAPP. 

The  NORSAR  and  NORESS/Finnmark  array  event  processors  differ  too  much 
to  have  a  common  program,  but  much  individual  coding  for  subprocesses 
is  common  to  all  packages.  The  main  topic  for  event  processor  develop¬ 
ment  will  be  code  structuring. 

Processes  involving  file  management,  plotting  and  computations  will  be 
separated  to  form  smaller  entitles.  This  is  done  both  in  order  to  get 
smaller  general  purpose  routines  and  to  identify  and  separate  codes 
chat  may  be  machine  dependent. 
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IV. 7 _ Hardware  Co  be  Installed  for  processing  of  data  from  the  new 

Finnmark,  northern  Norway,  regional  array 

Fig.  IV. 7.1  shows  how  the  data  processing  equipment  proposed  in  con¬ 
junction  with  the  newly  installed  regional  array  in  Finnmark,  Norway, 
will  be  Integrated  into  and  necworked  with  already  existing  resources 
at  NORSAR's  Data  Processing  Center  at  Kjeller.  The  data  acquisition 
systems  for  the  NORSAR  and  NORESS  arrays  (not  shown  in  this  figure) 
are  integrated  into  the  total  system  in  similar  ways. 

The  data  acquisition  and  real-time  processing  systems  (termed  FRS 
online  in  the  figure)  for  the  new  Finnmark  array  are  based  on  three 
Sun  3  computers.  The  bottom  Sun  3/260  is  interfaced  to  a  CIM 
(Communication  Interface  Module)  which  receives  data  from  the  Finnmark 
field  system  using  a  dedicated  satellite  link.  The  communication  pro¬ 
tocol  used  is  SDLC,  which  is  also  used  for  the  NORESS  installation. 

The  CIM  system  can  buffer  up  to  20  minutes  of  data.  In  normal  opera¬ 
tion,  data  are  sent  on  to  the  Sun  system  without  delays. 

Recording  of  both  high  frequency  (HF)  and  array  data  from  the  new 
Finnmark  site  will  utilize  disk  loops  that  can  hold  84  hours  of  array 
data  and  also  84  hours  of  HF  data.  One  disk  drive  will  be  assigned  to 
each  of  the  two  data  types.  The  Sun  3/260  computer  that  handles  the 
data  acquisition  tc*o-  will  also  be  running  a  state-of-health  program 
to  verify  data  quality  and  generate  statistics  on  uptime,  com¬ 
munication  errors  and  calibration  results.  All  data  will  be  per¬ 
manently  archived  on  tapes,  in  the  NORESS  format. 

Data  contained  in  the  disk  loops  will  be  generally  accessible  using 
NFS-  The  data  will  also  be  accessible  from  the  IBM  system  using 
NORSAR's  network  software  for  reading  files  that  are  stored  on  com¬ 
puters  in  the  network. 
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The  other  Sun  3/260  system  in  the  "FRS  online"-box  of  Fig.  IV. 7.1  will 
be  used  for  real-time  processing  and  subsequent  plotting  of  data. 
Relevant  plots  will  be  produced  on  tie  Versatec  system.  This  Sun  3/260 
computer  is  configured  in  exactly  the  same  way  as  the  acquisition  com¬ 
puter,  so  it  can  act  as  <x  backup  for  the  acquisition  task. 

The  Sun  3/110  will  be  usea  for  quality  checking  and  array  monitoring 
purposes . 

The  IBM/PC  is  the  system  currently  being  used  tor  controlling  critical 
tasks  in  the  N0RSAR/N0RESS  system.  This  system  automatically  calls  the 
operator  on  duty  in  case  of  system  failures.  The  FRS  online  system 
will  be  Integrated  into  this  alarm  system.  PC/NFS  will  be  used  to  read 
current  status  directly  from  files  stored  in  the  Sun  online  system. 

The  third  SUN  3/260  system  wiLl  be  used  for  data  exchange,  data  analy¬ 
sis  and  program  development. 

In  addition  to  these  four  Sun  computers,  our  hardware  proposal  con¬ 
tains  another  two  Sun  3/60  systems  and  a  Sun  3/50  computer  (only  one 
of  the  Sun  3/60  systems  is  shown  ,n  the  figure).  These  systems  will  be 
available  in  the  network  as  general  purpose  work  stations. 


J .  Fyen 
T.  Kveerna 
B.  Paulsen 

R.  Paulsen 

S.  Mykkeltveit 
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V.  MAINTENANCE  ACTIVITIES 

v . i  Activities  In  the  field  and  at  the  Maintenance  Center 

Much  time  has  been  spent  on  preparatory  work  In  connection  with  the 
new  array  In  Karasjok,  Finnmark,  and  finally  weeks  were  spent  on  the 
Installation  Itself. 

In  spite  of  this  extensive  task,  corrective  maintenance  has,  1  a  few 
cases,  been  carried  out  in  the  NORSAR  array.  Otherwise,  visits  have 
been  limited  to  occurrences  of  more  acute  character,  such  as  power 
failure,  localization  of  cables  which  otherwise  most  certainly  would 
have  been  damag  1  by  excavators. 

NORESS  has  been  visited  sad  corrective  actions  carried  out. 

A  cable  connection  between  the  HF-equlpment  and  the  multiplexer 
(Megamux)  has  been  replaced  by  a  fiber  optic  link. 

Bias jd  has  been  visited  and  corre:tlve  mainte  mce  carried  out. 

At  the  NDPC,  NORSAR  communication  circuits  have  been  regularly 
checked. 

The  NORSAR  subarray  electronics  have  been  checked  and  monitored 
regularly. 

A  summary  of  the  activities  performed  during  the  period  is  given  in 
Table  V.l.  For  completeness  we  list  activities  carried  out  for  all 
systems,  i.e.,  NORSAR,  NORESS,  Blasjd  and  the  new  array  in  Karas jok. 


O.A.  Hansen 
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Subarray/ 

area _ Task _  _ Date 

NORSAR  P.W.  Larsen  attended  a  project  meeting  14  Apr 

with  Sandla  in  connection  with  the  new 
Finnmark  array. 

Karasjok  P.W.  Larsen  participated  in  a  survey  in  15,16  Apr 

Karasjok  together  with  representatives 
from  Sandia. 

NORSAR  P.W.  Larsen  attended  a  meeting  at  NORSAR  28  Apr 

together  with  the  project  manager  and 
Barlindhaug  A/S,  the  consultants 

Oslo  NMC  staff  participated  in  a  fiber  optic  27-30  Apr 

course  - 

NORESS  Repaired  a  DHL  70  card  tor  site  AO  21  Apr 

NDPC  Daily  routine  check  of  SP/LP  data.  Weekly  April 

calibration  of  SP/LP  seismometers.  Daily 
check  of  communication  systems. 

04C  P-W.  Larsen  met  representatives  from  NT A/  2  May 

Hamar  and  road  administration  2  May  in 
the  04C  area  to  clarify  responsibilities 
in  connection  with  04C  communication  cable 
rerouting . 

06C  Replacement  of  the  'ertical  long  period  19,20,21 

centering  device  (ROD) ,  Incl.  adjustments  May 

NORESS  NMC  staff  assisted  representatives  from  19,20,21 

Sandia  Nat.  Lab.  (Dick  Kromer)  in  connec¬ 
tion  with  replacing  the  cable  between  HF- 
equipment  and  the  Megamux  with  a  fiber  optic 
link  connection. 

NMC  Continued  planning  and  preparation  in  con-  May 

nection  with  the  Finnmark  array. 

NDPC  Daily  check  of  SP/LP  dfta.  Weekly  calibra-  May 

tion  of  SP/LP  seismome! ers .  Daily  checks 
of  communication  systens. 
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Subarray/ 

area 

Task 

Date 

01B  (05) 

Cable  localized  in  connection  with  exca¬ 
vation  damage. 

16  June 

0 IB  (03) 

In  connection  with  trenching,  cable  pointed 
out . 

24  June 

02B  (04) 

Cable  pointed  cut. 

17  June 

04C  (02) 

Cable  pointed  out. 

19  June 

06C  (03) 

Cable  work. 

17,18,22 

June 

Karas jok 

Survey  In  connection  with  the  new  array. 

11-15  June 

Bias  j o 

Various  repair  and  maintenance  work  carried 
out . 

16-20  June 

NORESS 

A  fiber  connector  r  paired  at  site  C2 

27  June 

NDPC 

SP/LP  daily  checks.  Calibration  of  SP/LP 
seismometers  weekly.  Daily  control  of 
communication  systems. 

June 

02B  (tel.) 

Corrective  maintenance  performed. 

21,23  July 

04C  (03) 

Cable  localized. 

29  July 

06C  (04) 

Found  damaged  cable. 

30  July 

Karas jok 

P.W.  Larsen  attended  a  meeting  in  Alta 
regarding  contractual  relations. 

6  July 

Karas jok 

P.W.  Larsen  visited  Karasjok  in  connection 
with  start  of  the  construction  work  (main 
building  and  cable  trenching). 

14-17  July 

FINE- A 

Visited  and  expanded  the  array  with  5  new 
C-points . 

29  July 

NDPC 

Daily  check  of  SP/LP  data.  Weekly  calibra¬ 
tion  of  SP/LP  seismometers.  Daily  control 
of  the  communication  systems. 

July 

Task 


Dace 


Subarray/ 

area 


06C 


Visit  in  connection  with  damaged  cable.  7  Aug 


0  IB 


Visit  in  connection  with  damaged  cable.  12  Aug 


NDPC  Daily  check  of  SP/LP  data.  LP  seismometers  Aug 

continuously  checked  with  regard  to  mass 
position  (MP) ,  and  free  period  (FP) .  Ad¬ 
justment  carried  out  when  outside  specifica¬ 
tions.  Communication  systems  checked  daily. 

Karasjok  NMC  staff  engaged  in  the  new  Finnmark  array  Sept 

throughout  the  period. 

NDPC  In  addition  to  ordinary  daily/weekly  checks/  Sept 

tests,  the  A/D  converters  of  all  the  data 
acquisition  equipment  (SLEMs)  have  been 
checked  with  respect  to  their  ability  to 
reproduce  correct  numbers. 


Table  V.l  Activities  in  the  field  and  at  the  NORSAR  Maintenance 

Center,  including  NDPC  activities  related  to  the  NORSAR 
array,  1  Apr  -  30  Sep  1987. 


V . 2 _ Improvements  and  modifications  (NORSAR  array) 

No  changes  since  last  reporting  period. 

V .  3 _ Array  s  tat  us 

As  of  30  Sept  87  the  following  channels  deviated  from  tolerances: 

01A  01  8  Hz  filter 

02 

04  30  dB  attenuation 


24  - 


01B  05 

08 

04C  01 

08 

06C  04 

05  Broadband  filter  installed. 
08 


O-A.  Hansen 


-  25  - 


VI.  DOCUMENTATION  DEVELOPED 

Bache,  T.C.,  S.R.  Bratt  and  H.  Bungum  (1986):  High-frequency 
attenuation  along  5  teleseismic  paths  from  Central  Asia, 
Geophys.  J.R.  astr.  Soc.,  85,  505-522. 

Loughran,  L.B.  (ed.):  Semiannual  Technical  Summary,  1  October  1986 
31  March  1987,  NORSAR  Scl.  Rep.  No.  2-86/87,  Kjeller,  Norway. 
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VII  SUMMARY  OF  TECHNICAL  REPORTS  /  PAPERS  PREPARED 

VII. 1  A  scattering  model  of  regional  Pn  wave  propagation 

The  Pn  phase  is  crucial  for  the  detection  and  location  of  regional 
events.  It  is  therefore  important  to  determine  its  characteristic 
properties.  The  properties  of  Pn  are  well  understood  in  one¬ 
dimensional  crust-mantle  models:  Pn  is  an  ordinary  head  wave  asso¬ 
ciated  with  the  Moho  in  models  of  uniform  plane  layers  and  Pn  can  be 
interpreted  as  a  sum  of  whispering  gallery  waves  forming  a  so-called 
interference  head  wave  in  models  of  spherical  layers  and/or  models 
having  a  positive  velocity  gradient  below  the  Moho  (Menke  and 
Richards,  1980).  The  phenomenon  of  the  whispering  gallery  can  explain 
the  relatively  long  duration  and  high  frequencies  of  Pn  ac  teleseismic 
distances.  However,  it  cannot  explain  the  often  observed  relatively 
large  amplitudes  in  the  later  part  of  the  signal.  Similar  charac¬ 
teristics  are  observed  at  regional  distances.  To  determine  the  cause 
of  these  characteristics  we  have  analyzed  in  some  detail  the  NORESS 
records  of  Pn  from  a  suite  of  6  mining  explosions  in  the  Blasji#  area. 
The  explosion  site  is  about  300  km  from  NORESS  in  an  azimuth  direction 
of  240' •  A  typical  record  section  of  Pn  from  these  explosions  is  shown 
in  Fig.  VII. 1.1. 

It  can  be  seen  that  the  first  arrival  is  relatively  weak.  In  fact, 
this  arrival  is  easily  missed  for  small  events,  whence  Pn  is  often 
associated  with  the  dominant  part  of  the  wave  train.  We  have  applied 
wlde-band  f-k  analysis  (Kvaerna  and  Doornbos,  1986)  to  the  Pn  wave 
trains  from  all  6  events.  Slowness  solutions  as  a  function  of  fre¬ 
quency  are  summarized  in  the  velocity-azimuth  spectrum  of  Fig. 

VII- 1.2.  The  solutions  show  a  pronounced  and  consistent  anomaly  within 
the  range  2-4  Hz;  this  is  also  the  range  where  the  signal  has  its 
maximum  energy  (Fig.  VII. 1.3).  From  an  analysis  of  slowness  as  a  func¬ 
tion  of  time  (Fig.  VII. 1.4  displays  the  results  for  one  of  the  events) 
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we  infer  chat  the  first  arrival  has  a  slowness  and  azimuth  consistent 
with  Pn  in  a  one-dimensional  crust-mantle  model,  and  the  anomaly  is 
related  to  wave  energy  being  delayed  by  about  0.5-0. 6  seconds.  The 
anomaly  cannot  be  generated  locally  near  the  surface  since  this  would 
significantly  perturb  the  wave  front;  the  wavenumber  solution  shows 
that  the  wave  front  is  in  fact  very  nearly  plane.  The  observational 
results  suggest  that  these  "Pn"  waves  are  actually  the  result  of  scat¬ 
tering  at  depth.  Ray  tracing  backward  in  the  direction  given  by  the 
measured  slowness  leads  to  the  result  that  the  0.5-0. 6  s  time  delay  is 
explained  by  a  scattering  source  within  the  depth  range  30-40  km.  The 
Moho  her<>  (s  usually  taken  to  be  a*-  about  35  km  denth.  Hence  the 
observational  results  are  consistent  with  scattering  by  topographic 
relief  of  the  Moho.  An  Interesting  geological  aspect  is  that  the 
inferred  location  of  the  proposed  topographic  feature,  70  km  from 
N0RESS  in  the  southwest  direction,  coincides  with  the  border  of  the 
Oslo  Graben. 

Scattering  will  of  course  affect  all  wave6  interacting  with  a  rough 
Moho  discontinuity.  However,  the  scattered  waves  usually  arrive  in 
the  coda  of  a  relatively  strong  primary  wave.  In  contrast,  the  first 
arriving  Pn  is  relatively  weatt  due  to  the  smali  coefficient  of  refrac¬ 
tion  through  the  Moho,  and  scattering  due  to  topography  of  the  boun¬ 
dary  may  dominate  the  wave  train. 

To  illustrate  these  concepts,  we  can  apply  a  recently  developed  method 
for  modelling  scattering  by  topographic  relief  (Doornbos,  1988).  The 
first  applications  of  this  method  concerned  the  relatively  complicated 
case  of  a  solid-liquid  boundary.  We  have  done  some  numerical  experi¬ 
ments  with  .  liquid  layer  over  a  solid  half  space,  letting  velocities 
and  densities  correspond  to  values  in  the  crust  and  mantle.  One  of  the 
results  is  given  in  Fig.  VII. 1.5.  Here  we  have  assumed  a  boundary 


topography  characterized  by  a  cori elation  length  of  6  km  and  an 
average  height  of  1  km,  and  the  wave  field  is  taken  monochromatic  at 
2  Hz.  The  figure  shows  energy  flux  of  P  transmitted  upward  through  the 
boundary,  as  a  function  of  slowness  of  incident  P.  Two  modes  of  scat¬ 
tering  are  shown:  (1)  The  specular  flux  E°  in  the  direction  defined  by 
the  plane  wave-plane  interface  concept.  The  specular  flux  through  a 
rough  interface  is  reduced  with  respect  to  the  flux  through  a  plane. 
(2)  The  diffuse  flux  E^C  due  to  multiple  scattering  in  all  (upward) 
directions.  E^C  does  not  exist  for  a  plane  interface.  The  figure 
illustrates  well  the  sharp  increase  of  the  ratio  E^^/E°  as  the 
slowness  approaches  the  critical  value  corresponding  to  Pn ,  thus  sup¬ 
porting  in  a  qualitative  way  the  scattering  model  for  propagation  of 
this  wave. 

It  should  be  noted  that  E°  and  Esc  are  not  observable  parameters.  What 
can  be  inferred  is  the  flux  at  a  receiver  location  on  the  surface.  To 
model  this  we  need  to  know  the  areal  extent  of  topography,  and  the 
decay  factor  accounting  for  attenuation  of  Pn  along  the  Moho.  Clearly 
these  parameters  as  well  as  those  describing  the  topography  needed  to 
compute  Esc  and  E°  are  presently  not  or  only  poorly  constrained,  and 
further  experiments  are  needed  to  establish  useful  constraints.  Our 
preliminary  results  nevertheless  suggest  that  the  scattering  model  of 
Pn  wave  propagation  is  viable,  and  consequently  that  a  careful 
calibration  is  needed  before  using  this  phase  for  event  location  and 
velocity  determination  purposes. 

D.J.  Doornbos,  Univ.  of  Oslo 

T.  Kvsrna 
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Fig.  VII.  1.1  Noress  record  of  the  Pn  phase  from  an  explosion  in  the 
Blasj0  area.  The  top  trace  is  the  single  channel  obser¬ 
vation,  and  the  lower  trace  is  the  steered  beam,  both 
bandpass  filtered  between  2.0  Hz  and  4.0  Hz.  Note  the 
relatively  weak  first  arrival  at  about  2.0  seconds  com¬ 
pared  to  the  stronger  signal  arriving  0. 5-0.6  seconds 
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SLOWNESS  ANALYSIS  -  BLASJ0  EVENTS 

MOVING  FREQUENCY  BANDS  -  PN  PHASE 
HAND:  0.8  I'Z  STEP:  0.4  HZ 


0  2  <  S  S  10 
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Fig.  VII. 1.2  Velocity  and  azimuth  .stimates  for  the  Pn  phases  as  a 

function  of  frequency  The  six  events  were  processed  by 
the  wlde-band  method  n  0.8  Hz  wide  frequency  windows 
with  0.4  Hz  steps.  Between  2  and  4  Hz  the  deviation  from 
theoretical  azimuth  is.  between  10  and  15  degrees,  and 
the  apparent  velocity  is  significantly  lower  than 
expected  for  the  Pn  plase. 
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Fig.  VII. 1.3  This  figures  shows  the  uncorrected  power  density  spectra 
of  both  the  Pn  phase  and  the  preceding  noise  for  a  typi¬ 
cal  Bias jp  event. 
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i  SLOWNESS  ANALYSIS  -  BLASJ0  EVENTS 

!  MOVINC  TIME  WINDOWS  -  PN  PHASE 

I  WINDOW:  0.8  SEC  STEP.  0.2  SEC 


Slowness  estimates  for  one  of  the  Blasj0  events  as  a 
function  of  time.  The  data  were  processed  by  the  wide¬ 
band  method  in  the  frrquency  range  2. 0-4.0  Hz.  The 
traced  data  were  prea’lgned  to  compensate  for  offset 
across  the  array. 


Fig.  VII. 1.4 


Slowness  (s/km) 


Fig.  VII- 1.5  Energy  flax  at  2  Hz  through  rough  solid-liquid  interface 
with  the  velocity-density  structure  Vp+/“  =  6. 8/8.1  km/s 
Vs+-/~  =  4.5/0  kir/s,  “  2. 9/3. 4  g/cm  .  The  boundary 

roughness  is  characterized  by  an  average  height  of  1  km 
and  a  correlation  length  of  6  km.  The  scattered  flux  is 

upward  through  the  boundary. - — -  :  specular 

directin' •,  :  integrated  diffuse  flux. 
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VIX.2 _ Preliminary  tests  for  surface  waves  In  2-D  structures 


Introduction 

Surface  waves  at  short  periods  (T  <  20  s),  and  especially  the  Lg 
crustal  surface  waves,  usually  exhibit  complex  wavetrains.  For 
example,  relative  amplitudes  of  the  Lg  to  body-wave  phases  observed  at 
the  NORSAR  array  were  shown  to  be  strongly  path-dependent  (Kennett  et 
al,  1985),  and  large  transversal  motion  was  observed  on  Lg-phases  pro¬ 
duced  by  explosive  sources  in  the  North  Sea  (Kennett  and  Mykkeltveit, 
1984). 

Lateral  heterogeneities  are  thought  to  be  responsible  for  this 
complexity,  but  due  to  the  lack  of  methods  which  are  available  to 
model  surface  wave  propagation  in  laterally  heterogeneous  structures, 
little  of  the  complex  wavetrains  has  yet  been  explained.  In  order  to 
get  some  further  insight  into  these  surface  wavetrains,  we  have  ini¬ 
tiated  a  research  effort  aiming  at  comparing  surface  wave  recordings 
at  the  NORSAR  and  NORESS  arrays  with  results  of  numerical  modelling, 
using  a  coupled  mode  scheme.  The  preliminary  results  presented  here  do 
not  intend  to  model  a  realistic  situation,  but  to  check  the  applicabi¬ 
lity  of  the  method,  and  especially  t)  analyze  the  influence  of  its  two 
main  parameters. 

Method 

We  give  here  an  outline  of  a  coupled  local  mode  method,  the  details  of 
which  can  be  found  in  Maupin  (1987).  It  is  appropriate  for  propagation 
across  2-D  structures,  the  angle  of  Incidence  of  the  waves  upon  the 
structure  being  possibly  nonperpendi  :ular  to  the  symmetry  direction  y 
of  the  structure.  The  wavefield  is  decomposed  into  a  laterally  varying 
sum  of  the  local  modes: 
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where  u  is  a  displacement-stress  vector  of  the  elastic  wavelield, 
ur  are  displacement-stress  vectors  of  the  local  modes,  and  (kr,p,0)  is 
the  local  wave  vector  of  the  mode  r  in  the  cartesian  coordinate  system 
(x,y,z),  at  the  frequency  w. 

The  lateral  heterogeneity  introduces  a  lateral  variation  of  the  wave- 
numbers  kr,  as  well  as  energy  transfers  between  modes  expressed  by 
lateral  variations  of  the  amplitude  coefficients  cr.  The  lateral 
variations  ot  these  two  quantitief  satisty  the  equations: 


dkr 

- Frr 

dx 

(2) 

dcr  _  1  X 

-  “  l  Fqr  -  *  exp(i  j  (k^-kt)  dC)  cq 

dx  q*r  (kT-kf)  0 


where  Fqr  is  an  expression  involving  the  local  mode  displacements  and 
tensions,  combined  with  the  elastic  coefficients  and  density  in  inter¬ 
face  terms,  and  with  their  lateral  derivatives  in  an  integral  over 
depth. 

The  system  (2)  is  transformed  to  yield  a  first-order  equation  in  x  for 
reflection  and  transmission  matrices,  which  is  solved  numerically  in  x 
using  a  fourth-order  Runge-Kutta  scheme. 
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The  method  Is  tested  on  the  continental  margin  model  designed  by 
Badal  and  Seron  ( L 98 7 )  (Fig.  VII. 2.1)  to  calculate  Love  wave 
transmission  by  finite  element  method.  We  use  an  integration  step  in 
x  of  5  km  and  a  variable  integratior  step  in  z  of  the  form  a(n_1)  dzg 
for  the  n-th  step,  with  dzQ  =  0.5  kr.  and  a  =  1.02.  In  the  following 
two  paragraphs,  we  investigate  the  influence  on  the  coupling  matrices 
of  the  precision  with  which  the  local  modes  are  calculated,  and  of  the 
number  of  modes  used  in  the  representation  of  the  wavefield. 


Influence  of  the  number  of  zones 

The  local  modes  ur  appear  free  of  li  teral  differentiation  in  the 
expression  Fqr  of  equation  (2).  In  order  to  facilitate  their  com¬ 
putation,  it  is  thus  possible  not  to  calculate  them  at  each  integra¬ 
tion  step  in  x,  but  to  take  them  as  constant  inside  zones  where  the 
total  lateral  variation  of  the  structure  is  small  enough  for  their 
shape  not  to  vary  significantly.  The  coarser  a  zoning  we  define  across 
the  laterally  heterogeneous  structure,  the  less  accurate  local  dis¬ 
placement  functions  we  use  in  the  calculation  of  the  expression  Fqr. 
which  in  turn  result  in  less  accurate  coupling  matrices  and  local 
wavenumbers . 

In  order  to  determine  how  fine  the  zoning  must  be,  we  divide  the  con¬ 
tinental  margin  of  Fig.  VII. 2.1  inti  2,  5  and  10  zones.  The  coupling 
equation  is  integrated  with  the  thr.  e  different  zonings,  using  the 
local  modes  at  a  period  of  20  s,  calculated  at  the  center  of  each  zone 
by  a  classical  Haskell-Dunkin  method.  The  resulting  local  phase  velo¬ 
cities  of  the  Rayleigh  wave  fundamental  mode  are  shown  on  Fig. 

VII. 2. 2a  as  a  function  of  x.  The  phase  velocities  calculated  at  the 
center  of  each  zone  by  the  Haskell-Dunkin  method,  as  well  as  the  local 
phase  velocities  calculated  in  each  of  the  smaller'  zones  using  the 
variational  method,  are  shown  on  Fig.  VII. 2. 2b  for  comparison. 
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The  variational  method  is  designed  to  calculate  tirst-order  variations 
of  the  phase  velocity  associated  with  velocity  and  density  variations 
as  a  function  of  depth,  preferably  for  flat  layer  boundaries.  Equation 
(2)  treats  separately  the  effect  01  phase  velocities  of  sloping  layer 
boundaries  and  velocity  variations  inside  the  layers.  It  gives 
obviously  more  accurate  results  thin  the  variational  method  in  this 
continental  margin  model,  where  the  slope  of  the  layer  boundaries  is 
important . 

The  accuracy  of  the  phase  velocity  improves  significantly  when  using  5 
zones  Instead  of  2,  but  remains  stable  between  5  and  10  zones.  The 
number  of  zones  required  to  fit  a  given  velocity  curve  can  be  foreseen 
from  the  linear  trend  of  the  integrated  phase  velocity  inside  a  zone. 
The  proper  zoning  for  a  given  model  and  frequency  depends  on  the  cur¬ 
vatures  of  the  different  mode  phase  velocity  variations.  The  discre¬ 
pancy  observed  between  the  integrated  phase  velocity  and  the  zonal 
phase  velocities  around  x  =  30  km  nay  be  related  to  a  systematic 
discrepancy  between  the  vertically  varying  velocity  profile  and  the 
set  of  homogeneous  layers  required  by  the  Haskell-Dunkin  method  to 
calculate  the  zonal  phase  velocities. 

The  zoning  must  be  designed  to  calculate  accurately  not  only  the  phase 
velocities,  but  also  the  transmission  and  reflection  matrices.  Fig. 
VII. 2. 3  shows  the  transmission  matrices  for  the  6  first  modes  of 
Rayleigh  and  Love  waves  at  right  angle  to  the  continental  margin  when 
the  modes  are  calculated  at  20  s  period  for  2,  5  and  10  zones.  The 
zoning  has  a  negligible  influence  on  the  transmission  matrices.  The 
reflection  matrices  are  zero  whatever  the  zoning  for  this  model. 

When  the  propagation  is  not  at  right  angle  to  the  structure,  Love  and 
Rayleigh  waves  are  coupled.  As  the  dispersion  curves  of  Love  and 
Rayleigh  modes  of  the  same  harmonic  rank  are  usually  very  close  and 
may  even  cross  each  other,  small  relative  errors  in  k  produce  large 
errors  in  (k<l-kr),  leading  to  a  1<  ss  of  accuracy  in  the  Rayleigh-Love 
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coupling  terras.  The  accuracy  of  the  phase  velocity  has  therefore  a 
stronger  influence  on  the  coupling  matrices  in  that  case. 


Influence  of  the  mode  cut-off 

To  be  complete,  the  wavefield  representation  (1)  should  include  all 
the  surface  wave  modes  and  body-wave  terras,  or  the  infinite  set  of 
modes  of  a  structure  equivalent  to  the  spherical  Earth.  For  practical 
reasons,  we  need  to  restrict  the  set  of  modes  to  those  which  have 
somewhere  along  the  structure  a  significant  energy  in  the  frequency 
range  we  are  studying.  Except  in  particular  cases  of  well-defined 
waveguides,  the  coupling  between  modes  usually  decreases  gradually 
with  harmonic  rank  difference  and  the  mode  cut-off  is  often  chosen 
somewhat  arbitrarily.  Therefore,  it  is  important  to  know  which  ele¬ 
ments  in  the  coupling  matrices  are  influenced  by  the  drop-off  of  the 
higher  modes. 

Transmission  matrices  for  Love  wave  crossing  the  continental  margin 
at  right  angle  were  calculated  for  the  6,  8  and  10  first  modes  at  a 
period  of  10  s  (Fig.  VII. 2. 4).  The  influence  of  the  mode  cut-off  can 
essentially  be  inferred  from  on  the  last  row  and  column  of  the  matri¬ 
ces.  It  is  clear  for  example  that  the  addition  of  modes  9  and  10  has 
no  significant  influence  on  the  transmission  matrix  of  the  first  7 
modes,  including  the  6th  mode  which  is  very  strongly  coupled  to  neigh¬ 
boring  modes. 

The  important  point  to  note  here  is  that  the  energy  of  the  higher 
modes  not  accounted  for  is  mainly  fed  into  the  coupling  terms  of  the 
higher  modes  of  the  representation,  and  not  systematically  into  the 
diagonal  terms.  Therefore  the  cut-off  does  not  produce  an  upper  bias 
of  the  self-transmission  terms. 


Conclusion 


We  have  shown  that  the  coupled-loc il  modes  method  can  be  applied  to 
model  wave  propagation  in  realistic  structures,  needing  a  reasonable 
amount  of  lateral  zones  and  local  nodes  to  achieve  a  satisfactory  pre¬ 
cision  in  transmission  and  reflection  coefficients.  The  method  has 
acceptable  computing  times  (typicaLly  10  minutes  of  CPU  time  on  the 
IBM  4381/P02  of  NORSAR  for  the  examples  displayed  in  Fi&s.  VII. 2. 3  and 
VII- 2. 4).  Next,  we  intend  to  apply  this  method  across  laterally 
heterogeneous  structures  around  Scandinavia,  like  the  North  Sea  Graben 
and  the  Tornquist-Teisseyre  lineament,  and  compare  the  results  of  the 
modelling  with  surface  wave  data  recorded  at  the  NORSAR  and  NORESS 
arrays. 


V.  Maupin,  Postdoctorate  Fellow 
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Fig.  VXI.2.1  Continental  margin  model  after  Badal  &  Ser6n  (1987). 


□  zonal  phase  velocity 

+. _ +  integrated  phase  velocity  using  10  zones 

x - X  integrated  phase  velocity  U'  ing  5  zones 

• - -  integrated  phase  velocity  using  2  zones 


Fig.  VII. 2. 2b 
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Fig-  VII. 2. 3  Transmission  matrices  in  X  for  the  6  first  Rayleigh  and 
Love  modes  across  a  continental  margin  model,  at  20  6  of 
period,  using  2,  5  and  10  zones.  The  matrix  elements  are 
amplitude  transmission  coefficients  for  modes  which 
carry  a  unit  energy  flux. 
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VII. 3 _ False  alarm  statistics  an  1  threshold  Jeterminat Ion 

for  regional  event  detect  ion 

In  order  to  Improve  the  event  detection  performance  during  automatic 
processing  of  regional  array  data,  the  detection  thresholds  for  the 
different  beams  need  to  be  closely  examined.  In  this  context,  several 
factors  should  be  considered: 

I)  We  wish  to  make  the  detectiin  process  as  sensitive  as  possible 
with  emphasis  on  small  events,  i.e.,  we  want  to  operate  on 

a  low  threshold. 

II)  We  need  to  specify  how  many  false  alarms  per  time  unit  can  be 
accepted;  e.g.  ,  operation  oc  the  NORSAR  array  has  shown  that 
false  alarm  rates  of  up  tc  1C  %  are  acceptable. 

III)  Too  many  false  alarms  will  cause  the  phase  association  algorithm 
to  produce  fictitious  events. 

lv)  The  number  of  false  alarms  increases  with  the  number  of  beams. 

In  an  extended  beam  set,  we  should  only  Include  the  beams  that 
improve  the  detectability  without  significantly  increasing 
the  number  of  false  alarms. 

In  this  study  we  have  attempted  to  determine  the  thresholds  for  a 
large  beam  set  from  false  alarm  considerations.  The  procedure  lias  been 
as  follows: 

1)  Based  on  previous  studies  (Kvserna  and  Mykkeltvelt,  1986)  and 
the  NORSAR  staff's  experience  with  processing  of  regional 
array  data  (RONAPP),  define  the  filter  bands  that  appear  most 
appropriate  for  regional  event  detection.  The  resulting  twelve 
filters  are  given  in  Table  VII. 3.1. 
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For  each  filter  band,  bind  the  best  array  sub-geometry  and  ine 
corresponding  steer1"?  delays  for  coherent  beamf orming .  For 
detection  of  P-waves  we  have  attempted  to  have  maximum  3  deci¬ 
bels  signal  loss  due  to  missteering.  For  coherent  beam  detec¬ 
tion  of  S-waves,  we  have  also  applied  the  3  dB  criterion.  The 
resulting  P-wave  beam  deployment  spanned  the  velocity  space 
from  6.0  km/s  to  infinity,  whereas  the  velocity  space  from 
about  3.8  km/s  to  5.3  km/s  was  spanned  by  the  S-beams .  This 
resulted  in  72  P-beams  and  132  S-beams. 

For  each  filter  band,  find  tie  best  array  sub-geometry  for 
incoherent  beamforming  of  the  vertical  channels.  The  sub- 
geometries  were  chosen  from  ue  criterion  that  the  noise  should 
be  uncorrelated  between  chanrjfs.  Due  to  the  wide  response  pat¬ 
tern  of  the  incoherent  beams  (Ringdal  et  al,  1975),  only  one 
incoherent  beam  (with  infinite  velocity)  was  formed  for  each 
filter  band. 

For  each  filter  band,  form  ai  incoherent  beam  from  the 
eight  horizontal  channels. 

Under  the  assumption  that  the  SNR  is  lognormally  distributed 
during  noise  conditions,  it  cm  be  shown  that  the  cumulative 
distribution  of  number  of  detections  versus  detection  threshold, 
approximately  follows  a  straight  line  when  both  are  plotted 
with  logarithmic  axes  (see  Appendix  1).  In  this  study,  we  have 
selected  a  time  interval  (1987/ :072. 04. 30.0  - 
1987/072:16.00.00.0),  and  run  all  beams  with  a  low  threshold. 
From  the  cumulative  distributions  we  have  extrapolated  the 
noise  slopes  down  to  a  common  level,  in  this  case  10  detec¬ 
tions,  and  found  the  correspt  iding  detection  thresholds. 
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In  addition  to  the  obtained  threshold  values,  a  description  of  the 
beam  deployment  is  given  in  Table  VII. 3. 2a  and  Table  VII. 3. 2b- 

Fig.  VII. 3.1  (a-1)  show  the  cumula  ive  distributions  ior  the  coherent 
P-beams  and  the  incoherent  beams  m  ide  from  the  vertical  channels. 

The  dashed  tangent  lines  are  the  inferred  noise  slopes  extrapolated 
down  to  the  10  detections  level.  At  this  threshold  we  thus  expect  to 
have  10  false  alarms  ("noise  detections")  in  each  filter  band  during 
the  examined  time  interval. 

From  Fig.  VII. 3.1  we  can  see  that  he  threshold  differences  between 
the  coherent  and  incoherent  beams  ary  between  4.5  and  7  decibels. 
From  this  observation  we  can  inter  that  for  detection  of  phases  where 
we  consistently  can  achieve  coherent  beam  gain  exceeding  these 
threshold  differences,  coherent  be  un  detection  is  superior  to  inco¬ 
herent  detection.  From  NORESS  experience  (Kvterna  and  Mykkeltveit, 
1986),  coherent  beamtorming  of  In  phases  will  typically  result  in  a 
SNR  gains  between  8  and  15  dB,  thus  coherent  beams  are  clearly 
superior  for  P-phase  defection.  For  secondary  phases  like  Lg  where  we 
may  not  achieve  coherent  beam  gain  exceeding  the  threshold  differen¬ 
ces,  incoherent  beamtorming  may  on; perform  the  coherent  procedure. 

Another  interesring  feature  seen  f 'om  Fig.  VII. 3. li)  and  Fig. 

VII. 3. Ij),  is  the  large  number  of  detections  exceeding  the  obtained 
thresholds  for  the  two  filters  8.0-16.0  Hz  and  10.0-16-0  Hz.  These 
filter  bands  have  been  subjected  to  additional  detection  analysis 
of  a  24  hours  time  interval,  where  we  found  that  the  majority  of 
the  detections  exceeding  the  obtai  ied  thresholds  occurred  during 
working  hours.  A  viable  hypothesi  ;  Is  therefore  that  most  of  these 
detections  are  dee  to  nearby  man-made  activity. 
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Fig.  VII. 3. 2  (a-d)  show  the  cumulative  distributions  for  the  inco¬ 
herent  beams  from  both  the  vertical  md  horizontal  channels.  An 
interesting  feature  is  that  the  inferred  thresholds  from  the  inco¬ 
herent  beams  based  on  8  horizontal  channels,  is  not  higher  than  the 
thresholds  inferred  from  the  incoherent  beams  based  on  22  vertical 
channels.  This  observation  suggests  that  incoherent  beamforming  on 
a  mixture  of  horizontal  channels  more  effectively  reduces  the 
noise  variance  than  using  vertical  c  lannels  alone. 

The  number  of  noise  detections  (10)  chosen  as  a  common  limit  for  each 
filter  band  in  this  study  is  of  course  rather  arbitrary,  and  has  been 
made  mainly  for  the  purpose  of  achieving  a  common  basis  for  com¬ 
parison.  In  practical  operation,  it  may  be  desirable  to  operate  at 
different  farse  alarm  rates,  and  the  corresponding  thresholds  may 
easily  be  inferred,  e.g.,  to  reduce  the  number  of  detections  in  the 
8-0-16.0  Hz  filter  band  by  a  factor  of  two  (Fig.  VII. 3. li),  the 
thresholds  have  to  be  raised  by  about  3  decibels. 

It  should  be  emphasized  that  the  false  alarms  discussed  in  this 
paper  are  due  tc  the  stochastic  nature  of  the  noise,  and 
that  intervals  with  different  noise  ;haracteristics  (Fycn,  1986) 
may  give  deviations  from  these  curves.  Noise  bursts  and  close 
events  which  may  be  considered  as  false  alarms  do  not  follow 
the  lognormal  distribution.  To  avoii  detection  of  these  signals, 
the  thresholds  would  have  to  be  raised  above  the  levels  obtained 
in  this  study.  Fig.  VII. 3. 3  shows  a  nistogram  of  RONAPP  detections 
with  apparent  velocities  less  than  3.0  km/s  for  a  period  of  86  days. 
These  signals  are  mostly  generated  by  local  activity  in  the 
northeastern  direction  from  NORESS.  These  detections  can  easily  be 
separated  from  the  'interesting'  loc  il  and  regional  events  on  the 
basis  of  their  apparent  velocities.  The  incoherent  beams  that  mostly 
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detect  these  local  events,  are  also  our  best  tool  for  detecting 
regional  Lg  phases-  So  raising  the  threshold  to  avoid  the  low  velocity 
signals,  would  also  imply  missing  a  number  of  regional  Lg  detections. 
This  illustrates  some  of  the  conflicting  interests  that  has  to  be 
taken  into  account  when  determining  thresholds  and  detection  beams. 

In  conclusion,  the  thresholds  obtained  from  the  extrapolated  noise 
slopes,  indicate  the  lowest  operational  level.  If  in  the  detection 
process  we  want  to  exclude  signals  generated  by  nearby  activity  or 
local  noise  bursts,  the  thresholds  should  be  raised.  But  as  outlined 
earlier,  many  other  factors  have  to  be  considered  for  a  final  beam 
deployment  and  the  corresponding  thresholds. 


T .  Kveerna 
S.  Kibsgaard 
F.  Ringdal 
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Appendix  1 


We  assume  that  the  quantity  STA/LTA  (short-term  to  long-term  average) 
computed  at  regular  intervals  on  a  seismic  trace  (beam  or  single  chan¬ 
nel)  follows  a  lognormal  distribution  during  noise  conditions  (Ringdal 
et  al,  1975). 

Setting  U  =  log  (STA/LTA)  we  thus  get  the  following  probability  den¬ 
sity  of  this  variable 
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where  <t>(z)  denotes  the  standard  Gaussian  densiLy. 

The  probability  of  U  exceeding  a  given  threshold  value  becomes 
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where  <J(z)  denotes  the  standard  Gaussian  distribution  iunction. 


With  a  logarithmic  frequency  axis,  we  obtain  trom  (2)  the  following 
slope  S(u)  of  the  distribution  function  at  a  given  threshold  u 


S(u)  = 
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For  values  of  t  >  3,  we  may  use  the  following  approximation,  with  a 
relative  error  of  less  than  10  per  cent: 


1  -  <E>(t)  “  -  *  <t>(t) 
t 


(4) 


Applying  this  approximation  to  (3),  and  also  noting  that  p  =<  0,  since 
STA/LTA  fluctuates  around  1  during  noise  conditions,  we  obtain 


u  u 

S(u)  - - -  for  —  >  3 

<5  0 


(5) 


In  practical  detector  operation,  we  are  only  interested  in  thresholds 

u 

corresponding  to  large  values  of  -  ,  thus  the  approximation  (5)  is 

o 

applicable.  For  example,  a  threshold  value  of  u  =  3.5  •  a  corresponds 
to  l  false  alarm  in  4000  independent  STA/LTA  computations,  whereas 
u  »  4a  gives  1  false  alarm  per  30,000  such  computations.  A  practical 
operating  threshold  is  likely  to  b  -■  somewhere  between  these  values; 
thus  the  slope  (5)  will  be  approximately  constant  around  the  threshold 
value,  for  a  given  a. 
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Nr. 

Prototype 

Type 

Low 

High 

Order 

BP01 

BU 

BP 

1.0 

3.0 

3 

BP02 

BU 

BP 

1.5 

3.5 

3 

BP03 

BU 

BP 

2.0 

4.0 

3 

BP04 

BU 

BP 

2.5 

4.5 

3 

BP05 

BU 

BP 

3.0 

5.0 

3 

BP06 

BU 

BP 

3.5 

5.5 

3 

BP07 

BU 

BP 

5.0 

7.0 

3 

BP08 

BU 

BP 

6.5 

8.5 

3 

BP09 

BU 

BP 

8.0 

16.0 

3 

BP10 

BU 

BP 

10.0 

16-0 

3 

BP1 1 

BU 

BP 

1.0 

2.0 

2 

BP12 

BU 

BP 

2.0 

3.0 

2 

Table  VII. 3.1  This  table  show  the  filters  applied  in  the  experiment. 

All  were  recursive  But terworth  bandpass  filters  of 
order  2  or  order  3. 


Coherent  P-Beams 


Filter 

Configuration 

#  beam: 

Test(dB) 

Threshold(dB)  #  Detections 

1.0- 

3.0  Hz 

A0  CD 

1 

6.0 

12.0 

28 

1.5- 

3.5  Hz 

A0  CD 

1 

6.0 

11.0 

40 

2.0- 

4.0  Hz 

A0  CD 

10 

6.0 

12.7 

51 

2.5- 

A. 5  Hz 

A0  BCD 

10 

6.0 

12.4 

50 

3.0- 

5.0  Hz 

A0  BC 

7 

6.0 

11.4 

58 

3.5- 

5.5  Hz 

A0  BC 

9 

6.0 

10.9 

64 

5.0- 

7.0  Hz 

A0  BC 

9 

6.0 

11.7 

64 

6.5- 

8-5  Hz 

A0  BC 

9 

6.0 

11.5 

70 

8.0- 

16.0  Hz 

A0AB 

8 

4.0 

8.4 

100 

10.0- 

16.0  Hz 

A0AB 

8 

4.0 

8.5 

114 

1.0- 

2.0  Hz 

A0  CD 

1 

6.0 

12.7 

29 

2.0- 

3.0  Hz 

A0  CD 

1 

6.0 

12.7 

28 

Coherent  S-Beams 


Filter  Configuration  ♦  beams  Test(dB)  Threshold(dB)  ♦  Detections 


1.0-  3.0  Hz 

A0  CD 

6 

6.0 

12.7 

40 

1.5-  3.5  Hz 

A0  CD 

12 

6-0 

12-9 

49 

2.0-  4.0  Hz 

A0  CD 

12 

6.0 

12.9 

58 

2.5-  4.5  Hz 

AO  BCD 

12 

6.0 

12.5 

66 

3.0-  5.0  Hz 

AO  BC 

12 

6.0 

12.3 

58 

3.5-  5.5  Hz 

AO  BC 

12 

6.0 

12.0 

61 

5.0-  7.0  Hz 

AO  BC 

12 

6.0 

12.7 

45 

6.5-  8.5  Hz 

AOABC 

12 

6.0 

11.8 

72 

8.0-16.0  Hz 

AOAB 

12 

6.0 

9.0 

93 

10.0-16-0  Hz 

AOAB 

12 

6.0 

9.0 

114 

1.0-  2.0  Hz 

AO  CD 

6 

b  •  0 

13.2 

33 

2.0-  3.0  Hz 

AO  CD 

12 

6.0 

13.7 

49 

Table  VII. 3. 2a  In  addition  to  the  inferred  thresholds  and  the  number 


of  detections  exceeding  these  thresholds,  a  description 
of  the  coherent  beam  deployment  for  detection  of  P  and 
S-phases  is  shown  in  thi9  table.  The  respective  coloums 
describe  the  following  parameters: 


Filter  : 

Configuration: 


#  beams  : 

Test  : 

Threshold  : 

#  Detections  : 


The  investigated  filter  band. 

The  array  sub-geomery  applied  in  the 
beamforming.  AO  means  AOZ,  A  means  A- 
ring,  etc. 

Number  of  coherent  beams  within  the 
filter  band. 

The  experimental  threshold  in  dB. 

The  threshold  in  dB  inferred  form  the 
slope  of  the  cumulative  distributions. 
Number  of  detections  exceeding  the 
deduced  thresholds. 


Incoherent  beams  on  vertical  channels 


Filter  Configuration  #  beams  Test(dB)  Threshold(dB)  #  Detections 


1.0-  3.0 

Hz 

A0 

CD 

1 

1.9 

6.4 

48 

1.5-  3.5 

Hz 

A0 

CD 

1 

1.9 

6.6 

61 

2.0-  4.0 

Hz 

A0 

CD 

1 

1-9 

6.9 

44 

2.5-  4.5 

Hz 

A0  BCD 

1 

1.9 

6.9 

42 

3.0-  5.0 

Hz 

A0  BCD 

1 

1.9 

7.2 

33 

3.5-  5.5 

Hz 

AO  BCD 

1 

1.9 

6.8 

36 

5-0-  7.0 

Hz 

A0  BCD 

1.2 

4.7 

72 

6.5-  8.5 

Hz 

AO  BCD 

1 

1.2 

4.5 

95 

8.0-16.0 

Hz 

AO  BCD 

1 

1.2 

4.0 

77 

10.0-16-0 

Hz 

AO  BCD 

1 

1.2 

4.0 

80 

i  rt  n  rs 

L»<J- 

Hz 

AO 

nr\ 

1/17 

1 

1.9 

6-8 

36 

2.0-  3-0 

Hz 

AO 

CD 

1 

1.9 

7.9 

45 

Incoherent  beams  on  horizontal  channels 


Filter  Configuration  #  beams  Test(dB)  Threshold(dB)  #  Detections 


1.0-  3.0 

Hz 

E-W  N-S 

1 

1.9 

6.8 

35 

1.5-  3.5 

Hz 

E-W  N-S 

1 

1.9 

7.0 

50 

2.0-  4.0 

Hz 

E-W  N-S 

1 

1.9 

7.4 

40 

2.5-  4.5 

Hz 

E-W  N-S 

1 

1.9 

6-8 

48 

3.0-  5.0 

Hz 

E-W  N-S 

1 

1.9 

6.1 

54 

3.5-  5.5 

Hz 

E-W  N-S 

1 

1.9 

5.3 

52 

5.0-  7.0 

Hz 

E-W  N-S 

1 

1.2 

5.1 

44 

6.5-  8.5 

Hz 

E-W  N-S 

1 

1.2 

4.8 

81 

8.0-16.0 

Hz 

E-W  N-S 

1 

1.2 

3.7 

81 

10.0-16.0 

Hz 

E-W  N-S 

1 

1.2 

3.7 

92 

1.0-  2.0 

Hz 

E-W  N-S 

1 

1.9 

7.3 

23 

2.0-  3.0 

Hz 

E-W  N-S 

1 

1.9 

8.6 

29 

Table  VII. 

•  3.2b 

In  this  table 

we  have 

shown 

the  same  parameters 

as 

Table  VII. 

■  3.2a 

i,  but  now  for 

the  incoherent  beams. 

E-W  means 

the 

channels 

ACE, 

C2E,  C4E  and 

C7E . 

N-S  means 

the 

channe  Is 

AON, 

C2N,  C4N  and 

C7N. 

Number  of  detections.  Humber  of  detections 


56  - 


10  -30  Hz 


3) 


(.5  -  3  5  Hz 


b) 


SNR  (dB) 


2.0  -  4.0  Hz 


2.5  -  4.S  Hz 


d) 


SNR  (<1B) 


Fig.  VII. 3.1  These  figures  (a-1)  show  the  cumulative  distributions 

for  the  coherent  P-beams  (dashed-dotted  curves,/  and  for 
the  incoherent  beams  made  from  the  vertical  channels 
(solid  lines)  in  all  filter  bands  investigated.  The 
dashed  lines  indicate  the  slopes  of  the  distribution  of 
'noise'  detections. 


3.0  -  SO  Hz  g,  3.5  -  5.5  Hz 


SNR  (dBJ 


SNR  (dB) 


0.0  4.0  8  0  12.0  18.0  20.0 


SNR  (dB) 


SNR  (dB) 


Fig-  VII. 3. 2  This  figure  (a-d)  show  the  cumulative  distributions  for 
the  incoherent  beams  from  the  vertical  channels  (solid 
lines)  and  from  the  horizontal  channels  (dotted  lines) 
for  four  different  filter  bands.  The  dashed  lines  indi 
cate  the  noise  slopes  inferred  from  the  'noise'  distri¬ 
butions  of  the  'horizontal'  beams. 


Fig.  VII. 3. 3  This  histogram  displays  the  azimuthal  distribution  of 
RONAPP  detections  during  the  86-day  period  indicated, 
with  apparent  velocities  less  than  3.0  km/s.  Within  the 
investigated  time  interval  there  were  totally  10000 
detections.  3140  of  these  detections  had  an  apparent 
velocity  less  than  3.0  km/s. 


6i 


VII • 4 _ Initial  results  trom  analysis  of  data  recorded  at  the  new 

regional  array  In  Flnnmark,  Norway 

During  the  summer  and  fall  of  1987,  a  regional  array  was  installed 
near  the  town  of  Karasjok  in  the  county  of  Finnmark  in  northern 
Norway.  The  new  Finnmark  array  was  designed  to  he  as  closely  as 
possible  a  copy  of  the  NORESS  regional  array,  which  was  established  in 
southern  Norway  in  1984.  The  geometries  of  the  two  arrays  arn  th“’"e- 
t ore  practically  Identical  (deviations  between  corresponding  sensor 
positions  are  of  the  order  of  a  few  tens  of  meters,  due  to  local 
terrain  conditions),  and  the  data  outputs  are  the  same  for  the  two 
arrays.  Fig.  VII. 4.1  shows  the  location  of  the  two  regional  arrays 
in  Norway,  and  also  the  location  of  the  FINESA  regional  array  in 
Finland,  which  was  described  by  Korhonen  et  al  (1987).  The  geometry  of 
NORESS  (and  for  most  practical  purposes,  also  the  geometry  of  the 
Finnmark  array)  is  shown  in  Fig.  VII. 4. 2. 

Data  from  the  new  Finnmark  array  hav  ■  been  transmitted  continuously 
via  satellite  to  the  NORSAR  data  processing  center  at  Kjeller  since 
November  1,  1987.  The  data  are  subjected  to  automatic  detection 
processing,  with  a  beam  deployment  identical  to  the  one  used  for 
NORESS. 


In  the  following,  we  report  on  some  findings  resulting  from  analysis 
of  data  from  the  new  Finnmark  array.  It  should  be  emphasized  at  the 
outset  that  the  available  data  cover  no  more  than  a  two-week  period, 
and  that  a  comprehensive  assessment  of  the  capabilities  of  the  new 
arxay  must  awa't  the  collection  of  data  covering  a  longer  time  span. 

Noise  spectra 

Fig.  VII. 4. 3a  shows  corrected  noise  ipectra  for  altogether  17  elements 
ot  the  new  Finnmark  array  (the  verti  :al  sensors  at  AO  and  the  C-  and 
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D-rings,  see  Fig.  VII. A. 2),  taken  tt  00.00  GMT  on  day  315.  For  com¬ 
parison,  NORESS  noise  spectra  tor  the  same  number  of  channels  and 
taken  at  the  same  time,  are  shown  in  Fig.  VII. A. 3b.  From  these 
figures,  it  is  seen  that  below  2.0  Hz,  the  Finnmark  array  experiences 
a  higher  noise  level  than  NORESS,  whereas  above  2.0  Hz  the  Finnmark 
site  is  clearly  the  quieter. 

The  high  noise  level  at  low  frequencies  at  the  Finnmark  site  has  been 
confirmed  by  other  data  and  is  typically  even  higher  than  shown  in 
Fig.  VII. A. 3a.  At  the  time  of  these  spectra,  the  nearby  coast  of  Finn¬ 
mark  experienced  a  wind  force  A,  w  ilch  is  moderate.  The  noise  level  at 
these  low  frequencies  is  generally  believed  to  be  governed  by  the 
passage  of  major  weather  fronts  ov  :r  the  open  ocean.  Therefore,  it 
should  not  be  unexpected  to  find  t  le  higher  microseismlc  noise  levels 
at  the  Finnmark  array,  since  this  array  is  located  closer  to  the  coast 
than  the  NORESS  a  ray. 

For  the  frequency  range  above  2  Hz,  the  noise  level  at  the  Finnmark 
site  appears  to  be  3-5  dB  below  that  of  NORESS.  A  possible  explanation 
here  is  that  the  noise  in  this  band  is  lower  at  the  northern  site 
because  of  a  lower  population  density  and  also  lower  level  of  traffic 
and  industrial  activities,  compared  to  the  NORESS  site.  These  noise 
levels  must  also  be  rated  as  low  relative  to  year-round  averages  for 
NORESS,  as  investigated  by  Fyen  (1987). 

Figs.  VII. A. Aa  and  Ab  each  show  ten  uncorrected  spectra,  taken  hourly 
between  00.00  GMT  and  10.00  GMT  of  day  315,  for  the  Finnmark  and 
NORESS  arrays,  respectively-  Each  single  spectrum  represents  an 
average  of  17  spectra  for  the  vertical  sensors  of  A0,  the  C-  and 
D-rings.  The  NORESS  spectra  show  the  well-established  (Fyen,  1986a, b; 
1987)  difference  between  night-time  and  day-time  noise  characteristics 
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(particularly  around  6  Hz).  The  Finnmark  data  are  generally  below 
those  of  the  NORESS  site  for  frequencies  above  2  Hz.  The  two  or  three 
curves  with  the  higher  noise  power  in  Fig.  VII. 4.4a  represent  cases  of 
high  frequency  noise  bursts  at  the  Finnmark  site  during  daytime.  These 
bursts  are  visually  confirmed  by  careful  inspection  of  the 
seismograms.  A  more  comprehensive  study  is  needed  to  clarify  the  ori¬ 
gin  of  this  noise.  There  is  so  far,  however,  no  indication  of  constant 
noise  sources  like  power  plants  and  sawmills- 

It  should  be  emphasized  again  that  in  order  to  establish  reliable 
estimates  of  ambient  noise  levels  at  the  new  array  site,  studies  like 
those  undertaken  by  Fyen  ( 1986a, b;  1987)  for  NORESS  are  needed.  The 
material  analyzed  so  far,  however,  indicates  that  in  the  range  of  pri¬ 
mary  interest  to  regional  seismic  verification  (i.e.,  above  2  Hz),  the 
noise  level  at  the  new  Finnmark  array  site  is  generally  somewhat  lower 
than  the  NORESS  noise  level. 


Noise  suppression  by  beamforming 

The  NORESS  array  has  proved  very  proficient  in  the  enhancement  by 
bearaforming  of  the  SNR,  yielding  gains  that  are  often  of  the  order  of 
or  even  in  excess  of  /N  (N  being  the  number  of  sensors  used  in  the 
beamforming).  It  has  been  shown  that  this  success  is  largely  due  to 
the  highly  effective  noise  suppression  that  can  be  obtained  by 
selecting  appropriate  sub-geometries  tor  the  various  signal 
frequencies.  As  a  first  check  on  the  new  array's  capabilities  in  this 
regard,  noise  suppression  curves  wer  •  computed  and  compared  with 
corresponding  results  from  NORESS. 


In  Fig.  VII. A. 5a,  Che  noise  suppression  in  the  frequency  range  0-20  Hz 
for  vertical  beamforming  (no  shifts  introduced)  is  shown  for  three 
one-minute  intervals  taken  hourly  at  00.00,  01.00  and  02.00  GMT  on  day 
315.  The  sub-geometry  used  is  that  of  A0,  the  B-  and  C-ring  instru¬ 
ments  (13  sensors).  The  dots  in  this  figure  represent  average  values 
for  NORESS  that  are  taken  from  Fyen  (1986c).  The  /N  level  is  at  about 
-11  dB,  and  the  general  impression  ..eft  from  this  figure  is  that  the 
new  array  is  as  effective  in  suppressing  noise  as  NORESS,  for  this 
sub-geometry.  For  another  sub-geometry,  comprising  the  sensors  of  A0, 
the  C-  and  D-rings,  corresponding  results  are  given  in  Fig.  VII. A. 5b. 
Again,  we  see  that  the  noise  suppression  capability  is  comparable  to 
or  maybe  even  better  than  the  average  performance  of  NORESS.  This 
strongly  suggests  that  the  spatial  characteristics  (e.g.,  correlation 
lengths  vs.  frequency)  of  the  noire  field  are  very  similar  to  those 
found  at  NORESS.  It  has  previously  been  established  (Korhonen  et  al, 
1987)  that  the  NORESS  and  FINESA  arrays  exhibit  strong  similarities  in 
this  regard. 

Analysis  of  data  from  two  regional  events  located  at  the  Finnmark 
array 

As  examples  of  regional  events  recorded  on  the  Finnmark  array,  we  pre¬ 
sent  the  records  for  two  presumed  raining  explosions  in  the  Kola  penin¬ 
sula  of  the  USSR. 

The  C-ring  seismograms  for  the  first  event  are  shown  in  Fig.  VII. A. 6. 
The  event  occurred  at  67.6'N,  34.0'E  (according  to  the  University  of 
Helsinki  bulletin),  at  an  epicentral  distance  of  A 08  km  and  an  azimuth 
of  117.9"  .  The  phases  Pn,  Pg,  Sn  and  Lg  can  be  clearly  identified  by 
visual  Inspection.  These  phases  were  subjected  to  wide-band  slowness 
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analysis,  with  results  given  in  Figs.  VIi.4.7  and  VII. 4. 8.  We  see  that 
the  phase  velocities  derived  are  in  the  expected  range  for  these  pha¬ 
ses  and  that  deviations  from  the  "true"  azimuth  are  within  6-7*. 

Data  for  the  second  event  are  shown  in  Fig.  VII. 4. 9.  This  event  is 
located  at  68.1'N,  33.2'E,  at  an  epicentral  distance  of  349  km  and  an 
azimuth  of  113.7" -  Besides  the  phases  identified  for  the  first  event, 
we  now  also  see  a  clearly  developed  Kg  phase.  It  is  of  particular 
Interest  to  note  the  difference  between  the  two  events  in  this  regard, 
particularly  since  they  are  separated  by  not  more  than  about  60  km. 

The  occurrence  of  Rg  waves  in  the  records  for  events  of  epicentral 
distances  of  the  order  of  350  km  also  sharply  contrasts  what  we  have 
found  at  NORESS,  where  Rg  waves  are  never  observed  beyond  100  km 
distance.  The  results  of  the  wide-band  slowness  analysis  of  the  phases 
Pn,  Pg,  Lg  and  Rg  for  this  event  are  shown  in  Figs.  VII. 4. 10  and 
Vil.4.1i.  Again  we  see  that  the  phas ;  velocities  are  reasonable,  and 
the  azimuths  deviate  by  not  more  than  5"  from  the  "true"  value. 


Regional  event  detection 

An  initial  study  has  been  made  comparing  the  regional  event  detection 
performance  of  the  two  arrays  in  Norway.  A  two-week  period  (Oct  31 
Nov  18.  1987)  was  selected  for  this  purpose,  and  analysis  of  RONAPP 
processing  results  for  the  two  arrays  was  conducted.  For  both  arrays, 
the  beam  deployments  and  thresholds  were  identical,  and  the  same  as 
those  used  for  the  past  two  years  in  regular  NORESS  operation. 

Fig.  VII. 4. 12  shows  a  map  displaying  all  regional  events  located  by 
NORESS  during  this  time  period,  whereas  Fig.  VII. 4. 13  gives  a  similar 
map  for  the  Finnmark  array.  We  recall  that  in  order  for  one  array  to 


locate  a  regional  event,  at  least  two  phases  (P  and  Lg)  from  that 
array  must  be  detected  and  associated. 

Comparing  these  two  figures  is  quite  instructive,  and  probably  gives  a 
reliable  Impression  of  what  can  be  expected  during  long-term 
operation.  The  actual  number  of  located  events  is  similar  for  the  two 
arrays  (NORESS  152,  Finnmark  117).  However,  there  is  almost  no  overlap 
of  the  two  populations;  in  fact  only  8  events  were  located  by  both 
arrays.  The  large  majority  of  located  events  are  within  500  km  of  tue 
respective  arrays,  and  represent  in  most  cases  presumed  local  explo¬ 
sions  of  low  magnitude  (Ml  <  1.5).  Sites  where  such  explosions  are 
clustered  can  be  easily  identified  on  the  plots-  It  is  noteworthy  in 
particular  that  the  Finnmark  array  detects  and  locates  a  large  number 
of  mining  explosions  in  the  Kola  peninsula. 

Fig.  VII. 4. 14  shows  a  map  of  all  events  of  estimated  Ml  >  2-0  located 
hv  at  least  one  array.  In  those  cases  when  both  arrays  located  the 
same  event,  the  location  by  the  closest  array  was  chosen.  Events  with 
at  least  one  confirming  phase  (p  or  Lg)  from  the  other  array  are 
encircled.  Details  pertaining  to  the  figure  are  given  in  Table 
VII. 4.1. 

Compared  to  the  previous  figures,  it  is  clear  that  relaxing  the  cri¬ 
terion  for  "common"  events  to  requiring  only  one  confirming  phase  from 
the  other  array  significantly  Increases  the  overlap  of  the 
populations.  The  majority  of  events  in  Fig.  VII. 4. 14  are  thus  detected 
by  both  arrays.  It  is  particularly  interesting  to  observe  the  good 
performance  for  the  event  cluster  near  65'N  40’E,  which  is  at  a  con¬ 
siderable  distance  from  both  arrays  (700  and  1500  km,  respectively). 
These  events  were  in  the  magnitude  range  2. 5-2. 7  and  the  locations 
have  been  Independently  confirmed  by  the  Finnish  network. 


In  almost  all  cases  of  confirming  detections  by  one  phase  only,  this 
corresponded  to  a  P  phase  (Table  VII-4.1).  This  result  is  somewhat 
surprising  in  view  of  earlier  P  and  Lg  detection  studies  for  NORESS 
and  may  not  be  representative.  Further  assessment  of  the  joint  detec¬ 
tion  and  location  potencral  ot  ehe  two  arrays  will  require  a  more 
extensive  data  base,  and  will  be  the  subject  of  further  study. 

S.  Mykkeltveit 

F.  Rlngdal 

J .  Fyen 

T.  KvaBrna 


References 


Fyen,  J.  (1986a):  NORESS  noise  spect -al  studies,  preliminary  report. 
Semiannual  Technical  Summary,  l  Oct  1985  -  31  Mar  1986,  NORSAR 
Sci.  Rep.  No.  2-85/86,  Kjeller,  Norway. 


Fyen,  J.  (1986b):  NORESS  noise  spectral  studies  -  noise  level  charac¬ 
teristics.  Semiannual  Technical  Summary,  1  Apr  1986  -  30  Sept  1986, 
NORSAR  Sci.  Rep.  No.  1-86/87,  Kjeller,  Norway. 


Fyen,  J.  (1986c):  NORESS  noise  spectral  studies  -  beam  suppression. 
Semiannual  Technical  Summary,  1  Apr  1986  -  30  Sept  1986,  NORSAR 
Sci.  Rep-  No.  1-86/87,  Kjeller,  Norway. 


Fyen,  J.  (1987):  NORESS  noise  spectral  studies.  Noise  level  character¬ 
istics.  Semiannual  Technical  Summary,  1  Oct  1986  -  31  Mar  1987, 
NORSAR  Set.  Rep.  No.  2-86/87,  Kjeller,  Norway. 


Korhonen,  H.  ,  S.  Pirhonen,  F.  Rlngdal,  S.  Mykkeltveit,  T.  Kveerna, 
P.W.  Larsen  and  R.  Paulsen  (1987):  The  Finesa  array  and 
preliminary  results  of  data  anal /sis.  Univ.  of  Helsinki,  Inst, 
of  Seismology,  Report  S-16. 


NORESS 


Finnmark 


Reference  Array 


Total  no.  of  events  located 
by  reference  array _ _ 


152 


117 


No-  of  events  Ml  >  2.0 


31 


No-  of  events  Ml  >  2.0 
located  by  the  refer¬ 
ence  array,  detected 
or  not  detected  by  the 
other  array. 


Both  P  and  Lg  |  d 


P  only 


14 


Lg  only 


1 


Not  detected  I  8 


42 


8 


12 


1 


21 


Table  VII. 4.1  Statistics  of  detected  and  located  regional  events 
for  the  two  arrays  in  Norway  during  a  two-week  test 
period . 


i 


LATITUDE  (DEG  N) 


Fig.  V1X 


LONGITUDE  (DEG  E) 


87/!  I /IS  If.  <2  SORSAR 


.4.1  The  figure  shows  the  network  of  three  regional  arrays  in 
Fennoscandia.  1:  The  NORESS  array  in  southern  Norway; 

2:  The  new  array  In  Finnmark,  northern  Norway;  and 
3:  The  F1NESA  array  in  Finland. 
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LEGEND: 

•  VERTICAL  SHORT  PERIOD 

O  3-COMPONENT  SHORT  PERIOD 

▲  3-COMPONENT  BROAD  BAND 

AND  3-COMPONENT  SHORT  PERIOD 


Fig.  VII. A. 2  The  geometry  ol  the  NORESS  array.  The  geometry  of  the 
new  Finnmark  array  comes  very  close  to  being  identical 
to  that  of  NORESS;  deviations  between  corresponding  ele¬ 
ment  positions  are  of  the  order  of  some  tens  of  meters. 
The  channel  assignments  (vertical  only  vs.  three- 
component,  short  period  vs.  broadband)  are  identical  for 
the  two  arrays.  The  short  period  instrument  at  the 
center  of  the  array  is  denoted  AO. 


F INMARK  NOISE 


Fig-  VII. 4. 3a  Noise  spectra  corrects  1  for  system  response  fo 
Finnmark  array  for  1?  vertical  channels  at  AO, 
and  D-rings.  The  spectra  are  base  don  one  minu 
at  00.00  GMT  on  day  313.  The  power  density  is 


NORESS  NOISE 


0.(  1.0  10.0  100 
I  i/2vfl  09:31:56  FREOLCNCY 

Fig.  VII. 4. 3b  Same  as  Fig.  VII. 4. 3a,  but  for  NORESS  data  tak. 
00.15  GMT  on  day  313. 


FISN1ARK  AVERAGE  NOISE  QO-IO 


Fig.  VII. 4a  Uncorrected  noise  spectra  for  the  Finnmark  array  for  ten 
one-minute  intervals  taken  hourly  between  00.00  and 
10.00  GMT  on  day  315.  Each  spectrum  represents  an 
average  of  the  17  vertical  sensors  of  A0,  the  C-  and 
D-rings . 

N0RESS  AVERAGE  NOISE  00-10 


Fig.  VII. 4. 4b  Same  as  Fig.  VII. 4. 4. i,  but  tor  NORESS  data  taken  hourly 
between  UO.UO  and  10.00  GMT  on  day  315. 
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FIWftARK  CRING  BEAM  SUPPRESSION 


Fig.  VII. 4. 5a  Finnmark  array  noise  suppression  by  beanuormrug  lor  the  sub¬ 
geometry  comprising  the  AO,  B-  and  C-ring  sensors.  To  produce 
these  curves,  a  vertical  beam  is  formed  and  the  spectrum  for 
this  beam  is  divided  by  the  average  of  the  single  sensor 
spectra.  The  three  curves  result  from  one  minute  of  data 
taken  hourly  at  00.00,  01.00  and  02.00  GMT  on  day  315.  The 
dots  represent  typical  NORESS  noise  suppression  values  for 
this  sub-geometry.  The  horizontal  line  at  -11.1  dB  represents 
/N  suppression  for  13  sensors. 

FIWHARK  TELEV  BEAM  SUPPRESSION 


Fig.  VII. 4. 5b  Same  as  Fig.  VII. 4. 5a,  but  for  the  sub-geometry  comprising 

the  sensors  of  AO,  the  C-  and  D-rings.  The  horizontal  line  at 
-12.3  dB  represents  /N  suppression  for  17  sensors. 
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Fig.  VII. A.  7 
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Wide -band  slowness  analysis 
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Fig.  VII. 4. 8  Wide-band  slowness  spectra  for  the  Sn  phase  (top)  and  Lg 
phase  (bottom)  for  the  event  In  Fig.  VII. 4. 6. 


31  TMO 


sums 


Fig.  VII-4.9 
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Fig.  VII. h. 10  Wid^-band  slowness  spectra  for  the  Pn  phase  (top)  and  Pg 
phase  (bottom)  for  the  event  In  Fig.  VII. 4. y. 
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Wide-band  slowness  analysis 
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Fig-  VII. 4.11  Wide-band  slowness  spe  .tra  for  the  Lg  phase  (top)  and  Rg 
phase  (bottom)  for  the  event  in  Fig.  VII. 4. 9. 
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LONGITUDE  (DEG  E) 


37/11/24  14  IS  RORSAR 


Fig.  VII. 4. 13  Regional  events  located  by  the  Finrmark  array  during  a 
two-week  test  period- 


LONGITUDE  (DEG  E) 


B7/tt/2S  I2<3  NORSAR 


VII. 14  Regional  events  of  Ml  >  2.0  located  by  at  least  one  ot 
t b»_  two  arrays  in  Norway  during  a  two-week  test  period. 
Events  with  at  least  one  confirming  phase  from  the  other 
array  are  encircled. 
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Lg  Measurements  at  NORSAR 


The  objective  of  this  study  is  to  investigate  the  potential  of  using 
NORSAR  recordea  P  coda  and  Lg  waves  to  obtain  stable  magnitude 
estimates  of  large  underground  explosions,  thereby  improving  the 
possibility  of  obtaining  reliable  yield  estimates  from  teleseismic 
recordings . 


The  NORSAR  array  (Bungum  et  al,  1971)  has  been  in  operation  since 
1970.  Originally  comprising  22  subarrays,  with  a  total  of  132  short 
period  vertical  seismometers,  the  array  was  reduced  to  7  subarrays  (42 
SPZ  seismometers)  in  1976.  In  order  to  obtain  consistency  over  time, 
this  analysis  has  been  restricted  to  data  from  these  7  subarrays.  With 
the  exception  of  occasional  subarray  or  Instrument  outages,  this  pro¬ 
vides  for  a  stable,  high  quality  recording  system,  making  possible 
reliable  comparison  of  seismic  events  recorded  during  the  entire  time 
period . 


The  general  characteristics  of  NORSAR  recorded  P  coda  and  Lg  waves  for 
Semipalatlnsk  explosions  have  earlier  been  studied  by  Ringdal  (1983) 
and  Baumgardt  (1985).  Ringdal  (1983)  found  by  comparing  with  ISC 
reported  mb  that  the  use  of  Lg  based  magnitudes  effectively  eliminated 
the  systematic  P-wave  magnitude  differences  observed  at  NORSAR  for 
explosions  from  the  two  main  test  si<es  within  the  Semipalatlnsk  area 
(Shagan  River  and  Degelen  Mountains)  Gupta  et  al  (1985)  studied 
NORSAR  P  and  P  coda  recordings  from  NTS  explosions,  and  found  that  P 
coda  measurements  provided  significantly  more  precise  yields  than 
those  based  on  the  initial  P.  They  also  obtained  promising  results  in 
applying  P  coda  measurements  to  NORSAR  recordings  from  Semipalatinsk 
explosions . 
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The  propagation  path  from  Semipalatinsk  to  NORSAR  is  indicated  in  Fig. 
VII. 5.1,  where  also  a  map  of  the  NORSAR  array  is  shown.  In  spite  of 
the  large  epicentral  distance  (4200  km),  the  Lg  phase  can  be  iden¬ 
tified  for  most  explosions  of  >  5.0.  However,  as  shown  in  the 
example  of  Fig.  VII. 5. 2,  the  Lg  amplitude  is  usually  only  slightly 
higher  than  that  of  the  preceding  P  coda,  and  this  feature  is  largely 
independent  of  event  size. 

The  data  base  for  this  study  consists  of  72  presumed  explosions  from 
the  Shagan  River  area,  USSR,  covering  the  time  period  1972  to  1987. 

All  such  events  of  ISC  mb  >  5.0  were  included,  except  for  a  few  cases 
of  NORSAR  system  outages.  The  event  set  is  listed  in  Table  VII.5.1, 
together  with  ISC  or  PDE  magnitude  estimates  (mb)  as  well  as  parame¬ 
ters  estimated  in  this  study.  For  most  of  the  events,  very  accurate 
location  estimates  have  been  computed  by  Marshall  et  al  (1985),  and 
additional  such  data  have  been  provided  by  P-  Marshall  (personal 
communication).  Otherwise,  NEIS  or  NORSAR  location  estimates  are  used. 
In  addition,  this  study  has  also  cade  use  of  maximum  likelihood  mj, 
magnitudes  estimated  at  Blacknest  from  ISC  data  (P.  Marshall,  personal 
communication)  and  mb  (Lg)  estimaies  published  by  Nuttli  (1986). 


Parameter  estimation  methods 

Lg  and  P  coda  magnitudes  were  estimated  from  NORSAR  data  for  each 
event  using  three  different  methods  (for  illustration,  see  Fig. 

VII. 5. 2).  All  three  methods  were  applied  to  filtered  NORSAR  SP  chan¬ 
nels,  using  a  0. 6-3.0  Hz  Butterworth  bandpass  filter.  In  earlier  stu¬ 
dies  (Ringdal,  1933,1,  tills  frequency  band  has  been  found  to  retain  the 
main  P  and  Lg  energy  while  suppressing  microseismic  noise. 

The  first  method  pertains  to  the  Lg  phase  and  comprises  direct 
measurement  'in  individual  traces  ..f  the  largest  cycle  with  a  period 
close  to  l  Hz,  within  a  group  velocity  window  of  3.25-3.70  km/s.  (The 
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average  Lg  group  velocity  is  close  to  3.5  km/s.)  The  amplitude  (zero 
to  peak)  and  the  corresponding  period  were  measured  and  converted  to 
ground  motion  A^j  (i'th  event,  j'th  channel).  The  amplitudes  were 
then  corrected  for  dispersion,  geometrical  spreading  and  anelastic 
attenuation  as  described  by  Nuttli  (1986),  using  parameter  values 
(Qo  «  700,  E,  =  0.4)  derived  by  Nuttli  for  the  nearby  station  KON.  This 
resulted  in  an  Lg  magnitude  estimate  MLA^j  for  each  analyzed  trace. 

Repeating  this  process  for  all  center  instruments  of  the  7  operational 
NORSAR  subarrays  (deleting  faulty  channels)  the  NORSAR  Lg  magnitude 
MLA  was  then  estimated  as: 


MLAi  "  W  WLAi  J  (1) 

where  N  (N  4  7)  is  the  numV.i  of  chan  le  Is. 

The  second  and  third  methods  both  cor  prise  automatic  RMS  measurements 
of  individual  channels  and  pertain  tc  P  coda  and  Lg  magnitudes, 
respectively.  For  the  P  coda  measurement,  a  time  window  of  30  seconds 
is  selected,  starting  20  seconds  after  P  onset-  The  LG  window  is  2 
minutes  long,  covering  the  group  velocity  interval  of  3-67  to  3.33 
km/s.  In  addition,  a  30  second  noise  window  immediately  before  P  onset 
was  selected  for  each  processed  channel.  The  estimation  procedure 
described  in  the  following  was  applied  for  each  event  to  all  42 
currently  operational  NORSAR  SP  channels,  except  that  faulty  channels 
were  deleted. 

We  denote  by  spjg  the  k'th  sample  of  the  filtered  trace  of  event  i, 
Instrument  j.  We  define  the  log  mean  square  Lij  by 
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Lij 


(2) 


where  K  is  the  number  of  samples  in  the  selected  time  window.  The 
quantity  is  defined  as 


Lt 


A  L‘J 


(3) 


where  M(<42)  is  the  number  of  operative  channels- 

We  denote  uj^',  and  hj*'  the  values  of  Lj  when  computed  over  the 
noise  window  (30  seconds,  K  =  600),  P  coda  window  (30  seconds,  K  “ 
600)  and  Lg  window  (120  seconds,  K.  »  2400),  respectively. 

The  N0RSAR  P  coda  PMS  magnitudes  (MCR^)  and  Lg  RP1S  magnitudes  (MLR^) 
are  then  estimated  by  correcting  or  noise  and  adjusting  to  standard 
scales  as  follows: 


MCR^  =  0.5  •  log[exp(L^^)  -  exp(Li^) J  +  Bq 


(4) 


MLR|  =  0.5  *  log [exp( Ljh)  _  etp(LiN)j  +  BL 


(5) 
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The  constants  Bp  and  are  determined  by  the  constraints: 


)(MCRi 

-  MBMLEp)  =  0 

(6) 

I(MLRi 

-  MBLGp)  «  0 

(7) 

where  MBLGi  is  the  Lg  magnitude  give:  by  Nuttli  (1986)  of  the  i ' th 
event,  MBMLEp  is  the  maximum  likelihood  estimate  discussed  earlier 
and  the  sums  are  taken  over  all  common  events. 

We  also  computed  uncorrected  P  coda  and  Lg  magnitudes  for  purpose  of 
comparison.  These  were  obtained  by  deleting  the  noise  term  in  (4)  and 
(5). 


Data  analysis 

The  three  estimation  methods  described  above  were  applied  to  the 
events  in  the  data  set  of  Table  VII. 5.1.  The  procedure  to  compensate 
for  background  noise  level  had  negligible  effect  on  the  RMS  P  coda 
magnitudes,  since  the  SNR  in  these  cases  was  consistently  high. 
However,  for  some  of  the  smaller  events  (mb(ISC)  <  5.6),  the  effect  on 
RMS  Lg  was  significant,  and  the  resulting  estimates  must  be  considered 
more  uncertain  than  for  the  large  explosions.  This  is  illustrated  in 
Figs.  VII. 5. 3  and  VII. 5. 4,  which  shov  ,  respectively,  the  uncorrected 
and  corrected  RMS  Lg  estimates  plotted  against  the  NORSAR  amplitude 
based  Lg  magnitudes. 

Table  VII. 5. 2  gives  a  list  of  standard  deviation  of  magnitude  dif¬ 
ferences  for  all  combinations  of  sources.  We  note  that  the  NORSAR  P 
coda  magnitude  correlate  well  with  ISC  maximum  likelihood  mp,  (o  = 
0.06),  as  is  also  seen  from  the  plot  in  Fig.  VII. 5. 5.  In  fact,  the 
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correlation  is  about  as  good  as  between  the  two  ISC-based  measure¬ 
ments  - 

It  has  earlier  been  noted  by  many  authors  that  P  coda  provides  much 

more  stable  mb  estimates  than  the  initial  P  phase,  which  is  subject  to 

strong  focussing  or  defocussing  along  its  ray  path.  Such  effects  are 
to  a  large  extent  averaged  out  in  the  P  coda  by  various  scattering 
mechanisms.  Nevertheless,  our  analysis  has  shown  that  P  coda  ampli¬ 
tudes  across  NORSAR  are  positively  correlated  with  initial  P  amplitude 
patterns-  Thus,  not  all  of  the  P  : ocussing/def ocussing  effects  are 
eliminated,  and,  by  reciprocity,  locussing  at  the  source  would  be 
expected  to  significantly  affect  I*  coda  magnitudes,  especially  over  a 
larger  eplcentral  aita  than  considered  in  this  paper. 

From  Table  VII. 5. 2  we  also  note  t  iat  the  various  Lg  based  measures  show 

surprisingly  little  correlation  w  th  each  other,  although  the  RMS 

noise  correction  explains  much  of  the  scatter  between  the  two  Lg  esti¬ 
mates  based  on  NORSAR  data.  Fig.  VII. 5.6  shows  NORSAR  RMS  Lg  versus 
Nuttli  (1986)  mf,  (Lg)  for  common  events-  Even  at  high  magnitudes, 
where  NORSAR  noise  corrections  aru  insignificant,  there  is  a  large 
scatter.  The  reasons  for  this  lac;  of  consistency  need  to  be  further 
investigated,  but  it  must  be  noted  that  the  NORSAR  measurements  have 
the  advantage  of  being  based  on  a  system  stable  over  time,  and  with 
high  quality  digital  recording. 

Comparing  the  Lg  based  measurements  to  the  ISC  maximum  likelihcud  mb, 
it  appears  that  the  NORSAR  RMS  Lg  shows  best  correlation  (o  =  0.092). 
Fig.  VII. 5. 7  compares  these  two  magnitude  measures,  whereas  Fig. 

VII. 5. 8  shows  NORSAR  P  coda  versus  NORSAR  RMS  Lg  magnitudes.  Comparing 
these  two  figures,  it  is  evident  that  Fig.  VII. 5. 8  has  a  slightly 
larger  scatter  (a  -  0.101)  and,  moreover,  apparently  could  be  split 
into  tvo  subpopulat tons  with  paraLlel  trends. 
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Figs.  VII. 5. 9  and  VII. 5. 10  show  that  this  is  no  coincidence.  On  Fig. 
VII. 5. 9,  the  NORSAR  P  coda  -  RMS  Lg  magnitude  differences  are  plotted 
as  a  function  of  event  location  (Marshall  et  al  (1985)  and  NEIS  date). 
The  P  coda  -  Lg  statistics  are  clearly  region  dependent  within  the 
Shagan  River  area;  in  particular  the  northeast  part  show  consistently 
low  P  coda  magnitudes,  whereas  the  southwest  part  shows  P  coda  magni¬ 
tudes  consistently  higher  than  those  based  on  Lg.  Fig.  VII. 5. 10  shows 
that  the  same  trend  is  apparent  when  comparing  ISC  maximum  likelihood 
mjj  to  NORSAR  RMS  Lg  infinitudes.  It  is  clear  that  application  of 
regional  corrections  would  improve  the  mutual  consistency  of  these 
magnitudes . 

Marshall  et  al  (1985)  found  that  explosions  in  the  northeast  and 
southwest  portions  of  the  Shagan  River  area  produce  distinctly  dif¬ 
ferent  P  waveforms  when  recorded  at  the  UK  seismological  array  sta¬ 
tions,  suggesting  that  the  Shagan  Ri\ er  site  can  be  subdivided  into 
two  test  areas  characterized  by  dif  f*  rent  geophysical  properties.  Our 
results  support  this  suggestion.  Our  data  show  an  average  bias  mb(P)  - 
mb(Lg)  of  -0.059  ±  0.014  for  NE  Shagan,  and  0.112  ±  0-009  for  SW 
Shagan.  Thus,  if  we  attempt  to  explain  this  anomaly  as  resulting  from 
the  systematic  differences  in  P  recordings  only,  we  obtain  a  relative 
mb(P)  bias  of  about  0.17  mb  units  between  the  two  areas.  On  the  other 
hand,  the  possibility  of  an  mb(I'g)  bias  cannot  be  ruled  out  either. 

Fig.  VII. 5. 11  Illustrates  the  relative  occurrence  of  large  Semipala- 
tinsk  explosions  after  1976  as  a  fun< tion  of  NORSAR  RMS  Lg  magnitudes. 
The  plot,  which  has  a  nominal  vertical  scaling,  has  been  obtained  by 
adding  Gaussian  density  functions  with  standard  deviations  of  0.015 
magnitude  units  centered  at  each  observed  magnitude  value-  The  three 
clear  peaks  on  this  plot  might  Indicate  different  yield  categories.  It 
is  interesting  to  note  that  the  highest  magnitude  group  includes 
explosions  trom  both  the  southwest  and  the  northeast  areas  of  Shagan 
River,  whereas  similar  plots  based  on  P-wave  magnitudes  would  place 
all  of  the  largest  events  in  the  southwest  area.  The  rightmost  peak  ot 
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the  diagram  (magnitude  6.06)  would  correspond  to  a  yield  between  150 
and  200  kilotons  TNT  if  Nuttli's  (1986)  formulae  are  used.  With  the 
inherent  uncertainty  in  the  absolute  calibration  level,  this  is 
clearly  not  inconsistent  with  the  150  kiloton  upper  limit  of  the 
Threshold  Test  Ban  Treaty. 

In  conclusion,  N0RSAR  based  P  coda  and  Lg  magnitudes  appear  to  provide 
magnitude  estimates  of  large  Semi  >alatinsk  explosions  comparable  in 
stability  to  those  of  a  world-wid  ;  network.  There  are  indications  that 
the  Lg  based  measurements  avoid  some  of  the  bias  inherent  in  P-based 
magnitudes.  It  would  furthermore  be  reasonable  to  expect  that  such 
measurements  based  on  network  averaging  would  provide  even  better 
stability.  In  conjunction  with  calibration  data  to  obtain  accurate 
absolute  reference  to  yield,  such  data  would  appear  to  provide  fo- 
very  reliable  monitoring  of  a  threshold  treaty. 


F.  Ringdal 
B.K.  Hokland 
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ORIGIN 

DATE 

ORIGIN 

TIME 

MB 

***  LC 
MLA 

(AMP) *** 

N  STD 

***  LG (RMS)*** 
MLR  N  STD 

***  PCODA  **** 
MCR  N  STD 

11/02/72 

01.26.57.7 

6.1 

6  .  132 

7 

0.040 

6.121 

42 

0.069 

6.207 

42 

0.058 

12/10/72 

04.27.07.7 

6.0 

6.157 

7 

0.066 

6  .  120 

42 

0.068 

6.013 

42 

0.079 

07/23/73 

01.22.57.8 

6.1 

6  . 185 

7 

0.082 

6  . 199 

41 

0.078 

6.286 

41 

0.060 

12/14/73 

07.46.57.0 

5.8 

5.951 

7 

0 . 087 

5.876 

42 

0.063 

5.826 

39 

0.074 

10/16/74 

06.32.57.5 

5.5 

5.455 

7 

0.086 

5.413 

42 

0.062 

5.596 

42 

0.065 

i  ■>  /27  n* 

ce . .  8 

5.6 

5  .  849 

7 

0.054 

5.711 

42 

0.060 

5.512 

42 

0.093 

Q4/21/15 

05.36.57.3 

5.6 

5.677 

7 

0.060 

5.550 

42 

0.054 

5.732 

42 

0.072 

10/29/75 

04.46.57.5 

5.8 

5.755 

7 

0.100 

5.632 

42 

0.062 

5.733 

42 

0.057 

12/25/75 

05.16.57.0 

5.7 

5.817 

7 

0.053 

5.804 

42 

0.067 

0.000 

0 

0.000 

04/21/76 

05.02.57.4 

5.3 

0.000 

0 

0.000 

0.000 

0 

0.000 

5.226 

42 

0.062 

06/09/76 

03.02.57.6 

5.3 

5.462 

7 

0.091 

5.203 

42 

0.056 

5. 126 

42 

0.093 

07/04/76 

02.56.57.7 

5.8 

5.831 

7 

0.036 

5 .813 

42 

0.063 

5.909 

42 

0.075 

08/28/76 

02.56.57.5 

5.8 

5.762 

7 

0.085 

5 .737 

41 

0.069 

5.761 

41 

0.063 

05/29/77 

02.56.57.8 

5.8 

5 . 7b8 

7 

0.066 

5.669 

39 

0.061 

5.796 

41 

0.057 

06/29/77 

03.06.58.0 

5.3 

5.274 

7 

0.021 

5.070 

39 

0.054 

5 . 279 

40 

0.072 

09/05/77 

03.02.57.8 

5.8 

5.953 

7 

0.075 

5.891 

39 

0.067 

5.789 

40 

0.072 

10/29/77 

03.07.02.9 

5.6 

5.889 

7 

0.083 

5 . 784 

39 

0.061 

5.768 

41 

0.070 

06/11/78 

02.56.57.7 

5.9 

5 . 824 

7 

0.066 

5.752 

39 

0.059 

5.918 

39 

0.072 

07/05/78 

02.46.57.3 

5.8 

5.823 

7 

0.086 

5.789 

39 

0.057 

5.877 

39 

0.066 

08/29/78 

02.37.06.5 

5.9 

6.005 

7 

0.060 

6 . 006 

39 

0 . 057 

5.914 

39 

0.090 

09/15/78 

02.36.57 .3 

6.C 

5.948 

7 

0.093 

5.906 

38 

0.059 

5.979 

41 

0.077 

11/04/78 

05.05.57.5 

5.6 

5.906 

7 

0.039 

5.691 

39 

0.061 

5.631 

41 

0.069 

11/29/78 

04.33.02.9 

6.0 

5.982 

7 

0.086 

5.969 

39 

0.067 

6.004 

41 

0.070 

06/23/79 

02.56.57.6 

6.2 

6.088 

4 

0.084 

6.065 

21 

0.070 

6  .  174 

23 

0.057 

07/07/79 

03.46.57.4 

5.8 

5.996 

7 

0.091 

5.967 

39 

0.074 

5.850 

41 

0.071 

08/04/79 

03.56.57.2 

6  .  1 

6.105 

7 

0.032 

6.098 

39 

0.062 

6  .  170 

41 

0.076 

10/28/79 

03.16.56.9 

6.0 

6.060 

6 

0.056 

6.033 

39 

0.063 

5.936 

41 

0.085 

12/02/79 

04.36.57.5 

6.0 

5.962 

5 

0.046 

5.903 

33 

0.067 

6.074 

35 

0.050 

06/29/80 

02.32.57.7 

5.7 

5 .722 

4 

0.074 

s.e-’c 

o 

0.004 

: .  7  2  ; 

2  Z 

0 .  C40 

09/14/80 

02.42.39.3 

6.2 

0.000 

0 

0.000 

0.000 

6 

0.000 

6.222 

36 

0.049 

10/12/80 

03.34.14.1 

5.9 

5.946 

6 

0.043 

5.912 

39 

0.061 

5.840 

41 

0  .  u69 

12/14/80 

03.47.06.6 

5.9 

5.887 

5 

0.075 

5.911 

39 

0.045 

5.965 

41 

0.070 

12/27/80 

04.09.08.2 

5.9 

5.948 

5 

0.040 

5.908 

39 

0.076 

5.949 

4  1 

0.069 

03/29/81 

04.03.50.0 

5.6 

5.719 

5 

0.052 

5.521 

39 

0.053 

5.519 

41 

0.U/0 

04/22/81 

01.17.11.4 

6  .  0 

5.945 

5 

0.060 

5.888 

33 

0.046 

5.991 

35 

0.065 

05/27/81 

03.58.12.1 

5.5 

5.576 

4 

0.020 

5.442 

27 

0.061 

5.320 

29 

0.098 

09/13/81 

02.17.18.2 

6  .  1 

6.115 

6 

0.081 

6.092 

34 

0.076 

6.129 

35 

0.071 

10/18/81 

03.57.02.6 

6.1 

5.951 

7 

0.072 

5.976 

39 

0 . 057 

6.069 

35 

0.065 

11/29/81 

03.35.08.7 

5.7 

5.789 

5 

0 . 078 

5.566 

33 

0.062 

5.693 

35 

0.071 

12/27/81 

03.43.14.1 

6.2 

6  .  046 

6 

0.040 

6.056 

39 

0.074 

6.183 

41 

0.065 

04/25/82 

03.23.05.4 

6.  i 

6.074 

6 

n  .072 

6.056 

39 

0 .077 

6.111 

39 

0.071 

07/04/82 

01.17.14.4 

6  .  1 

0.000 

0 

0.000 

0.000 

0 

0.000 

6.078 

32 

0 . 089 

08/31/82 

01.31.00.5 

5.3 

0.000 

0 

0.000 

0 . 000 

0 

0 . 000 

5 .356 

39 

0 . 065 

12/05/82 

03.37.12.6 

6  .  1 

5  .989 

6 

0.088 

5.968 

37 

0 . 078 

6 . 076 

36 

0.060 

12/26/82 

C3. 35. 14.1 

5.7 

5.834 

5 

0.065 

5  .  653 

39 

0.061 

5.731 

29 

0.061 

06/12/83 

02.36.43.5 

6  .  1 

6.051 

5 

0 . 077 

6.070 

24 

0.071 

6.132 

24 

0.072 

10/06/83 

01.47.06.5 

6.0 

5.917 

4 

0.021 

5.876 

18 

0.054 

6.034 

18 

0.068 

10/26/83 

0’ .55.04.8 

C  1 

5 . 9?c 

f, 

0 .  h?s 

5 .  Sn3 

32 

0.055 

6.097 

32 

0 .  CCS 

11/20/83 

03.27.04.4 

5.5 

0.000 

0 

o.ooo 

0.000 

0 

0.000 

5.481 

32 

0.071 

Table  VI  1.5. 1  List  of  events  used  in  this  study.  The  table  incLudes 
date,  origin  time,  ISC/PDE)  as  well  as  NORSAR  based 
measurements  discussed  in  the  text  and  standard 
deviations  across  NORSAR.  Zero  entries  indicare  either 
no  NORSAR  data  •»vaiI.ApU*  or  SNR  too  low  tor  reliable 
.noHs>uroim*nt  • 
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ORIGIN 

DATE 

ORIGIN 

TIME 

MB 

***  LG ( AMP ) *** 
MLA  N  STD 

***  LG  ( RMS )  *  *  * 
MLR  N  STD 

***  PCODA  **** 
MCR  N  STD 

02/19/84 

03.57.03.4 

5.9 

5.764 

j 

0.053 

5.724 

27 

0.051 

5.889 

27 

0.057 

03/07/84 

02.39.06.4 

5.7 

5.775 

5 

0.075 

5.696 

27 

0.058 

5.632 

27 

0.073 

03/29/84 

05.19.08.2 

5.9 

5.905 

5 

0.035 

5.897 

27 

0.057 

5.889 

27 

0.071 

04/25/84 

01.09.03.5 

6.0 

5.905 

6 

0 . 072 

5.862 

33 

0.066 

5.972 

33 

0.071 

05/26/84 

03.13.12.4 

6.0 

6.075 

6 

0.063 

6.064 

33 

0.064 

5.965 

33 

0.115 

07/14/84 

01.09.10.5 

6.2 

6.089 

6 

0.073 

6  .  043 

33 

0.074 

6.148 

31 

0.088 

10/27/84 

01.50.10.6 

6.2 

6.101 

6 

0.077 

6.075 

32 

0.067 

6.254 

32 

0.073 

12/02/84 

03.19.06.3 

5.8 

5.973 

5 

0.141 

5.864 

27 

0.076 

5.747 

25 

0.081 

12/16/84 

03.55.02.7 

6 . 1 

6.055 

5 

0.093 

6.030 

27 

0.066 

6  .  113 

27 

0.068 

12/28/84 

03.50. 10.7 

6.0 

«  OQQ 

6 

0.057 

5.978 

33 

0.073 

6  .  147 

33 

0.073 

02/10/85 

03.27.07.6 

5.9 

5.898 

7 

0.070 

5.791 

38 

0.069 

5.947 

39 

0.068 

04/25/85 

00.57.06.5 

5.9 

5.979 

5 

0.078 

5.850 

27 

0.065 

5.936 

29 

0.069 

06/15/85 

00.57.00.7 

6.0 

6.024 

5 

0.075 

5.966 

28 

0.062 

6.099 

28 

0.069 

06/30/85 

02.39.02.7 

6.0 

5.983 

4 

0.068 

5.917 

28 

0.059 

6.058 

28 

0.065 

07/20/85 

00.53.14.5 

5.9 

5.884 

7 

0.076 

5.855 

36 

0.071 

5.912 

36 

0.068 

03/12/87 

01.57.17.0 

5.4 

5.385 

6 

0.081 

5.223 

31 

0.059 

5.437 

33 

0.068 

04/03/87 

01.17.08.0 

6.2 

6.062 

6 

0.079 

6.048 

31 

0.070 

6.245 

32 

0.052 

04/17/87 

01.03.04.0 

6.0 

5.928 

6 

0.072 

5.895 

31 

0.069 

6.077 

32 

0.060 

06/20/87 

00.53.04.0 

6.1 

5.992 

7 

0.117 

5.962 

34 

0.070 

6.124 

37 

0.092 

08/02/87 

00.58.08.0 

5.9 

5.945 

6 

0.103 

5.866 

30 

0.074 

5.885 

30 

0.074 

11/15/87 

03. 31. OB. 0 

6.0 

6.024 

7 

0.088 

5.959 

34 

0.063 

6.098 

36 

0.058 

12/13/87 

03.21.08.0 

6.0 

6.094 

6 

0.088 

6.066 

29 

0.080 

6.081 

30 

0.066 

12/27/87 

03.05.08.0 

6.0 

6 . 109 

6 

0.106 

6.032 

30 

0.083 

6.211 

30 

0.052 

AVERAGE  STANDARD  DEVIATIONS: 

0.070 

0.064 

0.069 

STANDARD 

DEVIATION  OF 

MEAN 

VALUES ; 

0.028 

0.011 

0.012 

Tabli'  VII. 5.  I  (coat.) 
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XSC/PDE 

ISC(M-L) 

MB(LG) 

PCODA 

LG ( RMS ) 

LG ( AMP ) 

XSC/PDE 

0.000 

0.066 

0.119 

0.068 

0.093 

0.109 

ISC(M-L) 

0.066 

0.000 

0.120 

0.060 

0.092 

0.136 

MB ( LG ) 

0.119 

0.  120 

0  000 

0.127 

0.  116 

0.  135 

PCODA 

0.068 

0.060 

0.  227 

0. 000 

0.101 

0.  130 

LG(RMS) 

0.093 

0.092 

0.116 

0.101 

0.000 

0.065 

LG ( AMP } 

0.109 

0.136 

0.135 

0.130 

0.065 

0.000 

Table  VII. 5-2  Standard  deviations  of  differences  between  various 
magnitude  estimate;  discussed  In  the  text- 


Kig.  VII. 5- l  Map  showing  the  propagation  path  from  Seralpalatinsk  to 
NORSAR  (distance  approximately  4200  km).  The  original 
NORSAR  array  conflguratt >n  la  also  displayed. 


Fig.  VII. 5. 5  Plot  of  NORSAR  P  co>la  magnitudes  versus  ISC  maximum 

likelihood  mb-  (Slope  restricted  to  1.00;  dotted  lines 
correspond  to  plus/minus  two  standard  deviations.) 


MAGNITUDE  COMPARISON 

Stia  jan  River  events 
SLOPE-  fOO  INTERC--0  01  STD~  0  116 


HORSAR  LG  RMS 


Fig.  VII. 5. 6  Plot  of  Nut t 1 1  (196.)  mb  (Lg)  versus  NOkSAR  Lg  RMS 

magnitudes.  The  slo >e  and  dotted  lines  are  defined  as 
in  Fig.  VII. 5. 5-  No  e  the  significant  scatter  even  at 
high  magnitudes. 


MACNITUL'E  COMPARISON 

Shagaii  River  events 
SLOPE  =  100  IN  TERC~-0. 05  STD=  0  092 


Fig.  VII. 5. 7  Plot  ot  ISC  maximum  Likelihood  mb  versus  NORSAR  Lg  RHS 

magnitudes.  The  slope  and  dotted  lines  are  defined  as  in 
Fig.  VII. 5. 5. 


MACNITUOE  COMPARISON 

Shag  in  River  events 
SLOPE =  t.00  ISTERC^-a  Oe  STD=  0  tot 


Fig.  VtL.S-d  Plot  of  NORSAR  RMS  P  c.da  mb  versus  NORSAR  RMS  Lg  magni¬ 
tudes.  The  slope  and  dotted  lines  are  defined  as  in 
Fig.  Vll.5.5. 


LATITUDE 

49.9  SO.O  50.1 


FREQUENCY 

0.0  0  2  0.4  0.6  0.8  1. 1 


MAGNITUDES 

Shagan  River  events 


NORSAR  LG  RMS 


VII-5.ll  ReLatlve  frequency  <>f  occurrence  of  large  Shag. in  River 
explosions  after  19^6,  based  on  estimated  NORSAR  Lg  RM 
magnitudes.  The  figure  has  been  generated  by  adding 
Gaussian  probability  density  functions  with  standard 
deviations  of  0.015  centered  at  each  magnitude  value. 
Vote  the  three  distinct  peaks  In  the  diagram. 


1 

i 


I 


-  103  - 


VII  .6 _ Towards  an  optimum  beam  deployment  for  NORESS; 

experiments  with  a  North  Sea  /  western  Norway  data  base 

Several  of  our  previous  Investigations  have  addressed  the  problem  of 
determining  the  number  and  type  of  beams  needed  for  adequate  detection 
of  signals  recorded  on  NORESS.  For  example,  Kvaerna  and  Mykkeltveit 
(1986)  arrived  at  a  recommendation  for  an  optimum  beam  deployment  for 
detection  of  regional  and  tcleseismic  P-waves- 

As  described  in  several  papers,  e.g.,  Mykkeltveit  and  Bungum  (1984), 
automatic  determinations  of  event  epicenters  using  our  on-line  pro¬ 
cessing  package  requires  detection  and  identification  of  both  P-  and 
S-arrivals  from  each  event.  This  study  utilizes  data  from  a  fairly 
extensive  data  base  of  regional  events  recorded  at  NORESS,  and  we 
attempt  to  determine  an  optimum  beam  deployment  with  respect  to  detec¬ 
tion  of  all  phases  from  these  events. 


Data  base  and  detection  processing 

The  data  base  comprises  101  events  tr >m  the  North  Sea  /  western  Norway 
i  region  with  locations  as  shown  in  Fig.  VII. 6.1.  These  events  were 

located  by  the  western  Norway  network  and/or  the  network  operated  by 
the  British  Geological  Survey.  The  eplcentral  distance  range  from 
NORESS  is  250  -  1100  km,  and  the  magnitudes  are  between  1.6  and  4.2. 

This  event  data  base  comprises  well  located  events  for  which  it  is 
possible  to  accurately  define  the  types  of  observed  phases  and  to 
read  their  onsets  manually,  as  illustrated  in  Fig.  VII. 6. 2. 


All  events  were  subjected  to  detection  processing  by  a  large  beam  set 
using  a  new  program  system  developed  at  N0RSAR  (Fyen,  1987).  Distri¬ 
buted  among  12  different  filter  bands,  this  beam  set  contained  72 
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steered  beams  for  detection  of  P-vaves,  132  beams  for  detection  of  S- 
waves ,  12  incoherent  beams  made  f>om  the  vertical  channels  and  12 
incoherent  beams  made  from  the  horizontal  channels.  In  selecting  this 
beam  set  we  attempted  to  include  all  candidates  that  could  conceivably 
contribute  to  increasing  the  detection  capability.  A  description  of 
this  beam  set  is  given  in  section  VII. 3  of  this  report  (Kvaerna  et  al, 
1987). 

The  output  from  the  detection  analysis  is  a  report  containing  parame¬ 
ters  such  as  an  identification  of  the  detecting  beam,  detection  time, 
STA,  LTA  and  STA/LTA  (SNR).  For  all  events,  the  detection  reports 
were  divided  into  Pn,  Pg,  Sn  and  i,g  detections.  Table  VII. 6.1  shows 
an  example  of  a  report  where  the  different  detections  are 
identified. 

The  essence  of  the  following  is  a  strategy  for  determination  of  a  beam 
set  that  meets  certain  performance  criteria  and  at  the  same  time  con¬ 
tains  the  lowest  possible  number  of  beams,  selected  from  the  set  of 
228  beams  described  above. 


Analysis  procedure 

In  the  work  presented  in  section  VII. 3  of  this  report  (Kveerna  et  al, 
1987),  we  attempted  to  infer  the  STA/LTA  (SNR)  detection  thresholds  ot 
the  respective  beams  from  noise  statistics.  For  each  phase  detected 
in  the  present  investigation,  these  thresholds  were  used  to  create  a 
basis  for  comparison  between  the  different  beams-  The  procedure  for 
comparing  the  maximum  SNR  values  can  be  illustrated  by  the  following 
example : 


i)  Assume  that  a  Pn  phase  Is  detected  by  three  different  beams, 

by  beam  A,  with  a  maximum  SNR  of  26.0  dB, 

by  beam  B,  with  a  maximum  SNR  of  23.5  dB  and 

by  beam  C,  with  a  maximum  SNR  of  22.0  dB. 

ii)  Assume  that  the  detection  thresholds  inferred  from  noise  statistics 
are , 

for  beam  A,  12.0  dB 

for  beam  B,  10.0  dB  and 

for  beam  C,  6-0  dB 

iii)  If  we  take  the  detection  thresholds  to  represent 

a  common  basis  for  comparison,  a  unified  quality  measure 
of  the  maximum  SNR  values  will  be  the  number  of  decibels 
exceeding  these  thresholds.  These  quality  measures,  from  now 
on  referred  to  as  'the  adjusted  SNR  values',  will  for  the 
three  beams  in  this  example  be; 

Adjusted  SNR  value  for  beam  A:  (26.0  -  12.0)  dB  =  14.0  dB. 

Adjusted  SNR  value  for  beam  B ;  (23.5  -  10.0)  dB  =  13-5  dB. 

Adjusted  SNR  value  for  beam  C:  (22.0  -  6.0)  dB  =  16-0  dB. 

From  this  we  conclude  that  beam  C  is  the  'best'  beam  for  detecting 
this  particular  Pn  phase. 

Rather  than  progressing  in  our  investigation  with  the  concept  of  a 
'best'  beam,  we  introduce  the  concept  of  a  'valid'  beam  or  'valid' 
detection:  The  maximum  obtainable  SNR  value  of  a  detection  will  be 
reduced  by  phenomena  like  mis-steering  of  the  coherent  beams,  random 
increase  in  the  noise  level  and  lowered  data  quality  on  single  chan¬ 
nels.  We  therefore  decided  that  all  ietections  with  adjusted  SNR 
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values  within  3  decibels  of  the  best  adjusted  SNK  for  the  phase,  were 
to  be  counted  as  'valid'  detections.  This  implies  that  each  detected 
phase  will  have  several  possible  candidates  for  the  'best'  beam. 

The  3  decibels  acceptance  level  i;;  chosen  rather  subjectively, 
however,  the  coherent  beam  deployment  is  sufficently  dense  that  signal 
losses  due  to  mis-steering  are  wi'hin  3  decibels.  In  later  analyses, 
this  acceptance  level  can  of  course  be  reduced  or  increased. 

Table  VII. 6. 2  shows  the  maximum  SNR  values  and  corresponding  adjusted 
SNR  values  for  one  of  the  Pn-phases  in  this  investigation,  for  all  48 
beam  types. 

Results 

After  finding  the  'valid'  detections  for  all  phases  in  the  data 
base  of  101  events,  we  were  in  a  position  to  evaluate  which  beams 
should  be  included  in  the  optimum  deployment.  Tables  VII .6 .( 3a- 3d) 
show  the  'valid'  detections  for  all  Pn,  Pg,  Sn  and  Lg-phases, 
respectively . 

When  considering  the  Pn  phase  in  Table  VII. 6. 3a,  we  found  that  beam 
number  10  had  the  maximum  number  of  'valid'  entries,  60  out  of  a  total 
number  of  98.  From  this  we  concluded  that  this  beam  should  obviously 
be  included  in  the  final  beam  deployment.  An  interesting  observation 
was  that  many  of  the  Pn  phases  from  the  North  Sea  /  western  Norway 
region  had  dominant  signal  frequencies  above  10  Hz. 

Beam  number  6  had  the  highest  nura  ier  of  'valid'  detections  for  the  Pg 
phase,  see  Table  VII. 6. 3b.  For  the  Sn  phase,  (Table  VII. 6. 3c),  beam 
number  45  or  46  was  an  obvious  candidate.  Beam  number  46  had  66  'va¬ 
lid'  entries,  which  is  about  80%  of  the  total  number  of  Sn  detections. 


Like  the  Pn  phase,  the  Sn  phases  from  this  region  contain  much  energy 
at  high  frequencies. 

For  the  Lg  phase  (Table  VII. 6. 3d)  no  specific  beam  is  an  obvious  can¬ 
didate  for  inclusion  in  the  recommended  beam  deployment.  Beam  42, 
however,  exhibits  the  highest  number  of  'valid'  entries. 

In  order  to  find  the  optimum  beam  deployment,  we  attempted  to  group 
the  beams  together  in  various  ways.  We  required  that  all  (or  almost 
all)  detected  phases  should  be  represented  by  at  least  one  'valid' 
beam  in  the  final  deployment.  First,  we  attempted  to  do  the  grouping 
for  each  phase  type  individually. 

For  the  Pn-phase,  the  number  of  ’valid'  entries  for  the  various 
beams  are  given  in  the  last  column  of  Table  VII. 6.3a.  We  first  select 
the  beam  with  the  highest  number  of  'valid'  entries,  in  this  case  beam 
number  10,  and  that  beam  handles  60  of  the  detected  Pn-phases.  Next, 
beam  number  6  accommodates  22  of  the  remaining  38  beams,  and  so  on. 

By  iteration,  we  ended  up  with  a  priority  list  of  beams  as  given  in 
Table  VII. 6. 4a. 

From  this  we  can  infer  that  8  beams  are  required  to  cover  all  Pn 
phases,  but  5  beams  alone  cover  96%  of  the  'valid'  entries.  Note 
that,  e.g.,  beam  10  (10.0  -  16.0  P)  includes  eight  coherent  beams  with 
different  steering  delays.  The  actual  azimuth  range  for  this  data  base 
is  180-360  degrees,  so  it  is  sufficient  to  distribute  about  4  coherent 
P-bearas  in  each  filter  band  to  span  this  area.  If  we  in  addition 
accept  that  it  is  sufficient  to  cover  96%  of  the  valid  detections,  the 
optimum  beam  deployment  for  detecting  Pn-phases  from  this  specific 
region  will  contain  16  coherent  P-beams  and  one  incoherent  beam  made 
from  the  vertical  channels.  If  we  require  the  coherent  beam  deployment 
to  span  all  azimuth  directions  (symmetrically),  about  32  coherent  P- 
beams  must  be  Included. 
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Tne  beam  priority  lists  for  the  P;>,  Sn  and  Lg  phases  are  given 
in  Tables  VII.6.4b-d,  respectively. 

An  interesting  observation  is  the  large  number  of  beams  required 
to  cover  the  Pg  detections.  From  Tables  VII.6.4a-d  we  can  infer  the 
beam  deployment  required  for  optimum  detection  under  the  given 
requirements,  when  treating  the  phases  individually;  see  Table 
VII. 6. 5. 


When  recommending  specific  beam  diployments,  it  is  important  to 
keep  in  mind  the  implications  on  •  omputational  load  for  the 
on-line  processor.  This  load  is  irimarily  determined  by  the 
total  number  of  traces  that  must  he  bandpass  filtered,  and  the 
number  of  beams  that  are  created. 

For  processing  the  data  with  the  deployment  of  Table  VII. 6. 5, 

178  single  channels  must  be  filtered,  and  55  beams  must  be  formed, 
each  of  which  is  subjected  to  STA/LTA-type  detection  processing. 

These  numbers  can  be  considerably  reduced  if  we  apply  an  algorithm 
where  all  phases  are  treated  at  the  same  time.  We  first  include 
the  most  obvious  beam  candidates  (beams  no.  46,  no.  10,  no.  42  and 
no.  6),  and  then  start  the  interative  sorting  procedure.  The  results 
from  this  procedure  are  given  in  Table  VII. 6. 6.  Compared  to  the  num¬ 
bers  in  Table  VII. 6. 5,  the  size  oi  the  beam  deployment  is  reduced 
by  30-402 . 

As  a  final  experiment,  we  decided  to  exclude  the  coherent  S-beams 
from  the  analysis.  The  results  are  given  in  Table  VII. 6. 7. 

Compared  to  the  results  given  in  Table  VII. 6. 6,  the  number  of  beams  to 
be  formed  were  marginally  reduced,  but  the  number  of  traces  to  be 
filtered  were  increased  due  to  including  more  incoherent  beams  made 
from  the  vertical  channels. 


-  109 


Conclusion 

For  detection  of  Pn-phases  from  the  North  Sea  /  western  Norway  region, 
the  high  frequency  coherent  P-beams  will  contribute  significantly  to 
an  optimized  beam  deployment.  For  the  more  low  frequency  Pg  phase, 
the  coherent  P-beam  filtered  between  3.5  and  5.5  Hz  seem  to  do  best. 

Including  incoherent  beams  made  from  the  horizontal  channels,  will 
improve  the  detectability  of  Sn  significantly.  We  also  found  that  the 
Sn  phase  in  most  cases  had  dominant  frequencies  around  10  Hz. 

Like  Pg,  the  Lg  phase  had  maximum  SNL  in  the  range  1.5-4. 5  Hz.  Inco¬ 
herent  beams  made  from  both  the  horizontal  and  vertical  channels  must 
be  included  for  optimum  detection  of  the  Lg  phase. 

A  tool  for  determining  the  optimum  beam  deployment  has  been  developed. 
To  validate  the  method  and  find  a  beam  deployment  that  adequately 
covers  also  other  relevant  regions,  we  plan  to  perform  similar  studies 
for  events  from  other  areas,  like  the  western  USSR. 

In  our  future  work,  then,  we  want  to  test  whether  the  symmetric  beam 
deployment  Inferred  from  the  North  Sea/western  Norway  data  base  is 
sufficient  for  optimum  detection  of  events  from  other  regions,  or 
if  other  beam  types  have  to  be  added  to  the  beam  deployment. 

T.  Kvasrna 
S.  Kibsgaard 
S.  Mykkeltveit 
F.  Ringdal 
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Manually  read  arrival  clava: 

Pn:  14.11.56.4  Pg:  14.12.05.5  Sn:  14:12.39.0 
Datactlon  results  for  steered  P-beaas. 


Tine  Interval  subjected  to  detection  processing: 


017:14.09.46.0  -  14.13.46.0 


Bean  Detection  interval 

STA 

LTA 

SNR  ^detection* 

***>  pn 

10906  017:14.11.56.6 

_ 

12.03.9 

827.01 

lb. 52 

44.661 

1 

H0807  017:14.11.56.9 

- 

12.06.4 

647.98 

12.91 

50.173 

4 

J1006  017:14.11.56.9 

- 

12.04.9 

543.86 

11.58 

46.969 

9 

E0505  017:14.11.57.1 

- 

12.03.4 

318.56 

26.33 

12.099 

2 

F0607  017:14.11.57.1 

- 

12.03.4 

454.06 

26.09 

17.406 

5 

C0707  017:14.11.57.1 

- 

12.03.6 

793.77 

30.41 

26.106 

5 

D0408  017:14.11.57.4 

- 

12.03.4 

263.52 

23.11 

11.403 

9 

CO 308  017:14.11.57.6 

- 

12.03.1 

279.83 

30.10 

9.297 

13 

B0 201  017:14.11.58.4 

- 

11.58. 9 

112.92 

27.44 

4.116 

1 

***>  Pg 

H0803  017:14.12.05.4 

_ 

12.11.9 

234.43 

41.59 

5.637 

16 

P0607  017:14.12.06.1 

- 

12.10.6 

812.29 

90.77 

8.949 

8 

E0505  017:14.12.06.1 

- 

12.09.1 

520.37 

76.49 

6.803 

15 

C0707  017:14.12.06.6 

- 

12.12.9 

749.01 

132.91 

5.635 

4 

J1002  017:14.12.07.6 

- 

12.08.1 

128.03 

30.39 

4.213 

8 

D0408  017:14.12.08.1 

- 

12.08.9 

345.48 

74.75 

4.622 

11 

F0602  017:14.12.13.6 

- 

12.13.9 

138.17 

34.16 

4.045 

1 

***>  ^ 

H0805  017:14.12.40.1 

_ 

12.42.1 

318.42 

77.59 

4.104 

2 

D0402  017:14.12.40.4 

- 

12.40.9 

231.42 

55.96 

4.136 

2 

G0702  017:14.12.40.4 

- 

12.41.4 

516.81 

98.40 

5.252 

3 

C0302  017:14.12.40.4 

- 

12.40.9 

253.58 

57.01 

4.448 

2 

***>  Lg 

F0603  017:14.12.51.9 

- 

12.54.1 

7*V,.74 

103.15 

7.065 

1 

E0502  017:14.12.52.1 

- 

12.54.1 

662.00 

78.63 

8.419 

12 

D0405  017:14.12.52.6 

- 

12- 54. 4 

464.31 

83.58 

5.555 

2 

C0302  017:14.12.53.1 

- 

12.59.6 

449.74 

78.16 

5.754 

3 

H0806  017:14.12.53.4 

- 

12.53.9 

474.39 

112.32 

4.223 

l 

F0602  017:14.12.55.9 

- 

12.56.6 

538.07 

95.86 

5.613 

2 

B0201  017:14.12.57.4 

- 

12.57.6 

305.59 

73.26 

4.171 

2 

D0402  017:14.12.58.1 

- 

12.58.6 

497.49 

116.50 

4.270 

i 

Table  VII-6.1  The  detection  report  shown  in  this  table  is  the  output 
from  detection  analysis  of  the  event  displayed  in  Fig. 
VII. 6. 2.  We  have  here  only  given  the  results  from  pro¬ 
cessing  the  3-minute  time  Interval  containing  the  Pn, 

Pg,  Sn  and  Lg  phases,  with  the  72  coherent  beams  with  P- 
wave  velocities.  The  beam  types  are  identified  from  the 
first  character  of  the  beam  name.  The  last  two  digits 
of  the  beam  name  identify  the  coherent  beam  with  maximum 
SNR  within  that  beam  type.  In  the  experiment  outlined 
in  this  paper,  we  give  results  for  each  coherent  beam 
type  rather  than  for  i  idividual  coherent  beams.  The  last 
column  of  the  detectlo  i  report  gives  the  number  of 
detections  for  the  bea  a  type  within  the  detection 
interval,  i.e.,  the  ti le  interval  for  which  the  SNR 
values  exceed  the  threthold. 
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Distance:  400.9  km. 

Azimuth 

:  229.0  degrees 

ML 

2.7 

B*  »m  no 

Bean  type 

SNR  In  d» 

Adjusted  SNR 

l 

1.0  -  3.0  P 

0.00 

/ 

0.00 

2 

1.5  -  3.5  p 

12.30 

/ 

1.29 

3 

2.0  -  4.0  F. 

19.37 

/ 

6.67 

4 

2.5  -  4.5  P : 

21.14 

/ 

8.74 

5 

3.0  -  5.0  P : 

21.66 

/ 

10.25 

6 

3.5  -  5-5  F. 

24.82 

/ 

13.91 

7 

5.0  -  7.0  P: 

28.34 

/ 

16.63 

8 

6.5  -  8.5  P: 

34.01 

/  * 

22.51 

9 

8.0  -  16.0  P: 

33.00 

/  * 

24.60 

10 

10.0  -  16  0  P: 

33.44 

z  . 

24.94 

11 

1.0  -  2.0  P: 

0.00 

/ 

0.00 

12 

2.0  -  3.0  P: 

0.00 

/ 

0.00 

13 

1.0  -  3.0  S: 

0.00 

/ 

0.00 

14 

1.5  -  3.5  S: 

0.00 

/ 

0.00 

15 

2.0  -  4.0  S: 

0.00 

/ 

0.00 

16 

2.5  -  4.5  S: 

0.00 

/ 

0.00 

17 

3.0  -  5.0  S: 

15.56 

/ 

3.26 

18 

3.5  -  5.5  S: 

19.41 

/ 

7.41 

19 

5.0  -  7.0  S: 

20.82 

/ 

8.12 

20 

6.5  -  8.5  £: 

26.52 

/ 

14.72 

21 

8.0  -  16.0  S: 

28.17 

/ 

19.17 

22 

10.0  -  16.0  J  : 

31.88 

/  * 

22. 8P 

23 

1.0  -  2.0  :  : 

0.00 

i 

0.00 

24 

2.0  -  3.0  S: 

0.00 

/ 

0.00 

25 

1.0  -  3.0  V: 

0.00 

/ 

0.00 

26 

1.5  -  3.5  V: 

7.24 

/ 

0.64 

27 

2.0  -  4.0  V: 

0.00 

/ 

U.00 

28 

2.5  -  4.5  V: 

0.00 

/ 

0.00 

29 

3.0  -  5.0  V: 

11.43 

/ 

4.22 

30 

3.5  -  5.5  V: 

14.12 

/ 

7.32 

31 

5.0  -  7.0  V: 

18.69 

/ 

13.99 

32 

6.5  -  8.5  V: 

23.17 

/ 

18.67 

33 

8.0  -  16.0  V: 

26.63 

/  * 

22.63 

34 

10.0  -  16.0  V: 

27.94 

/  * 

23.94 

35 

1.0  -  2.0  V: 

0.00 

/ 

0.00 

36 

2.0  -  3.0  V: 

7.60 

/ 

0.00 

37 

1.0  -  3.0  H: 

0.00 

/ 

0.00 

38 

1.5  -  3.5  H: 

0.00 

/ 

0.00 

39 

2.0  -  4.0  H*. 

0.00 

/ 

0.00 

40 

2.5  -  4.5  H: 

0.00 

/ 

0.00 

41 

3.0  -  5.0  H: 

0.00 

/ 

0.00 

42 

3.5  -  5.5  Hr 

17.13 

/ 

11.84 

43 

5.0  -  7.0  Hr 

15.46 

/ 

10.36 

44 

6.5  -  8-5  Hr 

19.34 

i 

14.54 

45 

8.0  -  16.0  H: 

20.00 

/ 

16.30 

46 

10. 0  -  16.0  H: 

20.82 

/ 

17.12 

47 

1.0  -  2.0  H: 

0.00 

/ 

0.00 

48 

2.0  -  3.0  Hr 

0.00 

/ 

0.00 

Table  VII. 6. 2  This  table  gives  the  maximum  SNR  values  and  the 

corresponding  adjusted  SNR  values  for  the  Pn  phase  shown 
in  Fig.  VII. 6.2.  The  beam  types  with  the  characters  'P' 
or  'S'  are  coherent  beams  with  P-  or  S-wave  velocities. 
The  character  'V'  indicates  incoherent  beams  made  from 
the  vertical  channels,  whereas  'H'  indicates  incoherent 
beams  made  from  the  horizontal  channels.  The  beam  types 
with  a  star  (  *  )  in  the  column  for  adjusted  SNR  are 
'valid'  beams,  i.e.,  within  3  decibels  of  the  beam  type 
of  highest  SNR.  Buam  types  with  the  value  0.0  have  not 
detected  the  phase 
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Table  VII. 6. 3  d)  This  table  shows  the  ’valid'  detections  for  all  Lg  phases  in  the  data  base. 
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a) 


Beam 

10 

10 . 0  -  16.0 

P 

60 

entries 

Beam 

6 

3-5  -  5.5 

P 

22 

entries 

Beam 

7 

5-0  -  7.0 

P 

5 

entries 

96Z 

Beam 

9 

8-0  -  16.0 

P 

4 

entries 

Beam 

22 

10-0  -  16.0 

V 

3 

entries 

Beam 

4 

2-5 

- 

4.5 

P 

2 

entries 

Beam 

5 

3.0 

- 

5-0 

P 

1 

entry 

4Z 

Beam 

27 

2-0 

- 

4.0 

V 

1 

ent  ry 

Total 

number  of 

detected 

Pn-phases 

98 

b) 

Beam 

6 

3.5 

- 

5.5 

P 

38 

entries 

Beam 

4 

2.5 

- 

4.5 

P 

8 

entries 

Beam 

16 

2.5 

- 

4.5 

s 

6 

entries 

Beam 

8 

6. 3 

- 

8.5 

P 

4 

entries 

Beam 

31 

5.0 

- 

7.0 

V 

4 

entries 

97X 

Beam 

13 

1.0 

- 

3.0 

s 

3 

entries 

Beam 

25 

1.5 

- 

3.5 

V 

3 

entries 

Beam 

18 

2.5 

- 

5.5 

s 

2 

entries 

Beam 

3 

2.0 

_ 

4.0 

p 

l 

entry 

3Z 

Beam 

20 

6.5 

8.5 

V 

1 

entry 

Total 

number  of 

detected 

Pg-phases 

70 

c) 

Beam 

46 

10.0 

- 

16.0 

H 

66 

entries 

Beam 

41 

3.5 

- 

5.5 

H 

7 

ent  rips 

Beam 

43 

5.0 

- 

7.0 

H 

3 

entries 

967. 

Beam 

32 

6.5 

- 

8.5 

V 

2 

entries 

Beam 

14 

1-5 

- 

3.5 

s 

1 

entry 

Beam 

19 

5.0 

- 

7.0 

s 

l 

entry 

Beam 

37 

1.0 

- 

3.0 

H 

l 

ent  ry 

Beam 

44 

6.5 

8.5 

H 

1 

entry 

Total 

number  of 

detected 

Si:~pbases 

82 

<0 

Beam 

42 

3-5 

- 

5.5 

H 

49 

ent  r ies 

Beam 

25 

1 .0 

- 

3.0 

V 

19 

ent  r le6 

Beam 

35 

1.0 

- 

2.0 

V 

7 

entries 

96Z 

Beam 

16 

2.5 

- 

4.5 

s 

6 

entries 

Beam 

27 

2.0 

“ 

4.0 

V 

3 

entries 

Beam 

20 

6.5 

- 

8.5 

s 

2 

entries 

4Z 

Beam 

37 

10.0 

“ 

16.0 

P 

\ 

entry 

Total  number  of  detected  L^-phaaes  :  87 


Table  VII. 6. 4  (a-d)  In  rhese  tab’"s  we  give  the  beam  priority  lists 
for  detection  of  Pn,  Pg,  Sn  and  Lg  phases  that  result 
from  applying  the  Iterative  procedure  to  each  phase 
type  individually. 


Number  of  traces 
to  be  filtered 


Number  of  beams 
to  be  formed 


Beam 

46 

LO.O 

- 

16.0 

H 

8 

1 

Beam 

10 

10.0 

- 

16.0 

P 

4 

4 

Beam 

42 

3-5 

- 

5-5 

H 

8 

1 

Beam 

6 

3.5 

- 

5-5 

P 

4 

4 

Beam 

41 

3.5 

- 

5.5 

H 

8 

1 

Beam 

35 

1.0 

- 

2.0 

V 

17 

1 

Beam 

16 

2.5 

- 

4.5 

S 

6 

0 

Beam 

7 

5-0 

- 

7.0 

P 

4 

4 

Beam 

r 

8-0 

- 

16.0 

P 

4 

4 

Beam 

8 

6.5 

- 

8.5 

P 

4 

4 

Beam 

31 

5.0 

- 

7.0 

V 

22 

1 

Beam 

22 

10.0 

- 

16.0 

V 

22 

1 

Beam 

13 

1.0 

3.0 

s 

6 

6 

Beam 

25 

1.5 

- 

3-5 

V 

17 

1 

Beam 

43 

5.0 

- 

7.0 

H 

8 

1 

Beam 

18 

2.5 

- 

5.5 

S 

6 

6 

Beam 

32 

6.5 

- 

8.5 

V 

22 

1 

Beam 

14 

1.5 

- 

3.5 

s 

8 

8 

178  55 

Symmetric  deployment  :  222  99 


Table  VII. 6. 5  From  the  results  produced  by  the  iterative  procedure 
with  each  phase  type  treated  individually,  the  beam 
deployment  given  in  this  table  was  inferred.  For  each 
phase  type,  this  beam  deployment  cover  96%  of  the  'va¬ 
lid'  detections.  We  have  also  given  the  number  of  tra¬ 
ces  to  be  filtered  and  the  number  of  beams  to  be  formed. 
To  cover  the  North  Sea  /  western  Norway  region  178  tra¬ 
ces  must  be  filtered  and  55  beams  formed.  To  span  all 
azimuth  directions  by  the  coherent  beam  deployment  the 
numbers  will  be  222  and  99,  respectively.  Note  that  for 
incoherent  beamforming,  each  individual  trace  must  be 
filtered  prior  to  tie  beamforraing. 
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Humber  of  traces  Number  of  beams 
to  be  filtered  to  be  formed 


Beam 

46 

10.0 

- 

16.0 

H 

8 

1 

Beam 

10 

10.0 

- 

16.0 

P 

4 

4 

Beam 

42 

3.5 

- 

5.5 

H 

8 

1 

Beam 

6 

3.5 

- 

5.5 

P 

4 

4 

Beam 

25 

1.0 

- 

3.0 

V 

17 

1 

Beam 

4 

2.5 

- 

4.5 

P 

4 

4 

Beam 

8 

6.5 

- 

8.5 

P 

4 

4 

Beam 

39 

2.0 

- 

4.0 

H 

8 

1 

Beam 

35 

1.0 

- 

2.0 

V 

17 

1 

Beam 

18 

3.5 

- 

5.5 

s 

6 

6 

Beam 

31 

5.0 

- 

7.0 

V 

22 

1 

Beam 

22 

10.0 

- 

16.0 

V 

22 

1 

Beam 

7 

5.0 

- 

7.0 

P 

4 

4 

128  33 

Symmetric  deployment  :  134  59 


Table  VI 1.6 -6  The  beam  deployment  of  this  table  was  produced  by  the 
iterative  procedure  when  treating  all  phases  at  the 
same  tLrae,  after  first  having  required  that  beam  nos. 
46,  10,  42  and  6  should  be  included  in  the  final 
deployment.  Compared  to  the  results  presented  in  Table 
VII. 6. 5,  the  volume  of  the  beam  deployment  is  con¬ 
siderably  reduced. 
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Number  of  traces  Number  of  beams 
to  be  filtered  to  be  formed 


Beam 

46 

10.0 

- 

16.0 

H 

8 

1 

Beam 

10 

10.0 

16.0 

P 

4 

4 

Beam 

42 

3.5 

- 

5-5 

H 

8 

1 

Beam 

6 

3.5 

- 

5.5 

P 

4 

4 

Beam 

25 

1.0 

- 

3.0 

V 

17 

1 

Beam 

4 

2.5 

- 

4.5 

P 

4 

4 

Beam 

8 

6.5 

- 

8.5 

P 

4 

H 

Beam 

39 

2.0 

- 

4.0 

H 

8 

1 

Beam 

35 

1.0 

- 

2.0 

V 

17 

1 

Beam 

31 

5.0 

- 

7.0 

V 

22 

1 

Beam 

22 

10.0 

- 

16.0 

V 

22 

1 

Beam 

7 

5.0 

- 

7.0 

P 

4 

4 

Beam 

9 

8.0 

- 

16.0 

P 

4 

4 

Beam 

29 

3.0 

- 

5.0 

V 

22 

1 

148  32 

Symmetric  deployment  :  172  5h 


Tabic  VII. 6-7  The  procedure  for  finding  the  results  given  in  this 
table,  is  essentially  the  same  as  outlined  for  Table 
VII. 6. 6-  The  only  difference  is  that  we  have  not  allowed 
any  S-beams  to  be  included  in  the  final  beam  deployment. 


Fig.  VII. 6.1  The  figure  shows  the  locations  of  the  101  events  in  the 
North  Sea  /  western  Norway  data  base.  The  loaction  of 
the  NORESS  array  is  indicated  by  a  filled  circle. 
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Fig.  VII. 6. 2  To  identify  the  different  regional  phases  and  pick 

their  arrival  times,  several  filters  were  applied  to 
both  the  vertical  and  horizontal  components.  This 
figure  shows  filtered  NORESS  recordings  from  an  event  in 
southwestern  Norway.  The  distance  trom  NORESS  is  about 
400  km  and  the  local  magnitude  is  2.7. 
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VII .7 _ Local  geology  of  the  regional  array  sites  In  Norway 

Information  on  local  geological  conditions  at  the  NORESS  and  Finnmark 
array  sites  is  available  from  various  sources.  In  this  contribution, 
we  have  attempted  to  compile  information  that  is  considered  relevant 
from  the  seismo logical  viewpoint.  In  particular,  such  information  may 
be  used  to  provide  tentative  explanations  of  certain  observed  features 
in  the  seismograms. 

Local  geology  of  the  NORESS  site 

The  NORESS  site  is  located  within  the  Precambrian  (also  referred  to 
as  pre-Eocambrian)  rock  complexes  of  southeastern  Norway.  Maps  showing 
the  local  geology  of  the  siting  arer  have  not  yet  been  published. 
However,  geologists  of  the  Norwegiai  Geological  Survey  have  conducted 
surveys  in  the  area  and  have  made  available  to  us  hand-colored  maps  of 
the  local  surface  geology.  These  mafs  show  that  three  different  rock 
types  are  present  at  the  array  site,  namely:  (1)  granite,  (2)  gabbro 
and  (3)  rhyolite.  The  information  in  these  maps  is  transcribed  onto 
the  array  configuration  map  of  Fig.  VII. 7.1. 

Johnson  and  Olheft  (1984)  and  Bott  (1982)  give  the  following  typical 
values  for  density  and  P-wave  velocities  (at  a  pressure  of  0.3  Kbar, 
corresponding  to  a  depth  of  1  km)  for  these  rock  types: 

Density  (g/cm3)  Velocity  (km/s) 

Granite  2.66  5.91 

Gabbro  2.99  6.74 

Rhyolite  2.(0 

One  should  keep  In  mind  that  these  values  represent  averages  over  a 
large  number  of  samples  from  various  geological  environments.  Wwe 
think,  however,  that  these  values  could  be  used,  with  due  care,  for 


the  site  considered  here.  For  insiance,  such  values  can  be  used  to 
estimate  reflection  coefficients  if  possible  subsurface  interfaces. 

The  possible  existence  of  such  interfaces  beneath  the  array  can  be 
inferred  from  data  from  a  reflection  seismic  survey  that  was  carried 
out  in  1983.  Reflection  data  were  collected  along  the  profiling  line 
shown  in  Fig-  VII. 7.1.  The  length  of  the  line  was  approximately  2500  m 
and  charges  of  1.4  kg  of  explosives  were  fired  in  8  m  deep  boreholes 
spaced  at  40  m  intervals  along  the  profile.  Each  shot  was  recorded  on 
a  spread  of  24  geophones  at  20  m  intervals,  with  a  constant  distance 
of  120  m  between  the  shot  and  the  nearest  geophone.  The  data  were 
recorded  using  NORSAR's  reflection  profiling  instrumentation  (Texas 
Instruments  DFS-V) .  A  sampling  rate  of  2  ms  was  chosen,  and  the  analog 
filters  had  cutoffs  at  8  and  256  Hz.  The  recording  interval  was  20  s 
for  each  shot.  The  field  work  and  subsequent  processing  of  the  data 
was  organized  and  carried  out  by  N0RSAR  personnel,  notably  Havar 
Gj0ystdal,  Paul  W.  Larsen  and  Ottar  A.  Sandvin. 

The  data  were  subjected  to  conventional  CDP-stacking  with  a  10  m 
interval  between  consecutive  traces  in  the  stacked  section,  giving  a 
6-fold  coverage  (except  at  the  ends  of  the  profile).  The  section  shown 
in  Fig.  VII. 7. 2  resulted  from  this  processing  and  shows  the  data 
filtered  in  the  100-250  Hz  band.  The  section  covers  a  two-way  travel 
time  (TWT )  of  0. 0-1.0  s. 

Several  interesting  features  can  he  observed  in  the  section  of  Fig. 
VII. 7. 2.  First  of  all,  we  see  a  dipping  interface  that  appears  to 
intersect  the  surface  near  the  southern  end  of  the  profile.  The  TWT  of 
0.37  s  observed  at  the  northern  end  of  the  profile  indicates  a  depth 
of  roughly  1  km  at  that  point,  giving  a  dip  of  approximately  20’. 

Based  on  the  available  Information  on  the  surface  geology,  we  specu¬ 
late  that  this  horizon  represents  the  interface  between  the  granite 
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and  the  gabbro,  which  is  outcropping  near  the  southern  end  of  the  pro¬ 
file,  see  Fig.  VII. 7. 2.  The  reflection  coefficient  resulting  from  den¬ 
sity  and  velocity  changes  across  a  grani te/gabbro  interface  could 
become  significant  and  would  be  of  the  order  of  0.12,  using  the  values 
given  above. 

In  addition  to  this  dipping  interface,  there  are  at  least  two  clearly 
observable  events  in  Fig.  VII. 7. 2.  One  of  these  is  dipping  across  the 
middle  of  the  section  at  am  .id  0.5-0. 6  s  TWT,  and  the  other  is  an 
apparently  flat  horizon  at  about  0.8  s  TWT.  Suffice  it  here  to  point 
to  these  arrivals,  without  going  any  further  into  geological  interpre¬ 
tations  . 

The  recording  interval  of  20  s  for  each  shot  made  it  possible  to  look 
for  coherent  arrivals  from  depths  even  larger  than  the  base  of  the 
crust.  Extensive  plotting  of  sections  for  TWT  intervals  within  the 
1.0-20.0  s  range  only  revealed  one  such  arrival  that  could  be  traced 
over  any  length  of  distance  and  tentatively  attributed  to  a  definite 
structure  at  depth.  This  event  is  seen  at  about  5.5  s  TWT  in  the 
middle  of  the  section  of  Fig.  VII. 7. 3.  The  corresponding  depth  would 
be  approximately  18  km. 

Local  geology  of  the  Finnmark  array  site 

Only  very  limited  information  is  yet  ivailable  on  the  local  geology  of 
the  new  regional  array  site  in  Finnmark,  northern  Norway.  The  Nor¬ 
wegian  Geological  Survey  has  published  a  regional  geological  map  of 
the  Karas jok,  Finnmark  area  (Skalvoll,  1972).  This  map  shows  that  all 
seismometers  of  the  new  array  are  within  an  outcropping  gabbro  for¬ 
mation  approximately  U  x  10  km  large.  This  gabbro  window  is  located 
within  the  Karas jok  Group,  which  consists  of  metamorphic  sedimentary 
and  igenous  rocks  of  age  2000  million  years.  It  is  therefore  reason¬ 
able  to  assume  the  existence  of  high-density,  high-velocity  material 
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at  the  surface  aLl  across  the  new  array.  The  possible  existence  of 
subsurface  layering  is,  however,  mknown  at  this  time. 


S.  Mykkeltveit 
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Fig.  VII. 7.1  Configuration  ar.d  local  geology  of  the  NORESS  array.  The 
four  three-component  stations  are  marked  with  special 
symbols.  The  station  at  the  center  of  the  array  is  named 
AO.  The  surface  geology  of  the  site  is  transcribed  from 
a  map  provided  by  the  Norwegian  Geological  Survey.  A  60 
m  deep  borehole  was  drilled  at  site  AO.  According  to  the 
drilling  log,  this  hole  was  drilled  in  gabbro  all  the 
way  down  to  its  bottom.  A  reflection  seismic  survey  was 
conducted  with  shots  a  \d  geophones  along  the  line 
running  into  the  array  from  the  north. 


Fig.  VII. 7. 2  CDP-stacked  section,  for  the  TWT  interval  0. 0-1.0  s.  The  filter  used  was  100-250  Hz. 
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